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ABSTRACT 

A well-known connection between n strings winding around a circle and permutations 
of n objects plays a fundamental role in the string theory of large N two dimensional Yang 
Mills theory and elsewhere in topological and physical string theories. Basic questions in 
the enumeration of Feynman graphs can be expressed elegantly in terms of permutation 
groups. We show that these permutation techniques for Feynman graph enumeration, 
along with the Burnside counting lemma, lead to equalities between counting problems of 
Feynman graphs in scalar field theories and Quantum Electrodynamics with the counting 
of amplitudes in a string theory with torus or cylinder target space. This string theory 
arises in the large N expansion of two dimensional Yang Mills and is closely related to 
lattice gauge theory with S n gauge group. We collect and extend results on generating 
functions for Feynman graph counting, which connect directly with the string picture. 
We propose that the connection between string combinatorics and permutations has im- 
plications for QFT-string dualities, beyond the framework of large N gauge theory. 
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1 Introduction 



There is a basic connection between winding states of strings and partitions of numbers. If 
we have a string wound around a circle, it can have a winding number which is a positive 
integer. For multiple strings, the winding numbers can be added to define a total winding 
number. In the sector of the multi-string Hilbert space with total winding number n, 
the partitions of n, denoted p(n), correspond to the different types of string states. The 
number p(n) also plays an important role in connection with the symmetric group S n , of 
order n\, consisting of the re- arrangements of the integers {1, 2 • • • , n}. It is the number 
of conjugacy classes in S n . 

In generic string theories, a string state is characterized by a winding number along 
with additional vibrational and momentum quantum numbers. In simple cases, the wind- 
ing numbers completely characterize the string state. A particularly striking example of 
this simplification occurs in the string dual of two dimensional Yang Mills theory (2dYM) 
on a Riemann surface Eg with genus G. This theory, for U (N) gauge group, can be solved 
exactly pQ. The coefficients in the large N expansion of amplitudes in 2dYM defined on 
Riemann surfaces can be expressed in terms of sums over symmetric group elements. This 
combinatoric data has an interpretation in terms of branched covers of by a string 
worldsheets S/j, Riemann surfaces of genus h related to the order in the 1/N expansion 
[21 H]. This provides a beautiful realization of the ideas of [5]. These spaces of covers 
are also called Hurwitz spaces and they are spaces of holomorphic maps from to S^. 
These spaces are reviewed from the point of view of 2dYM in [6] . 

Of particular interest are manifolds £g,b, with genus G and B boundaries. The 
observables of 2dYM are defined as functions of the boundary holonomies taking values 
in the gauge group U(N). The gauge- invariant functions of these group elements are 
multi-traces, which are also classified for fixed degree n, by partitions. This leads to 
an interpretation of the boundary partition functions in terms of covering spaces of the 
manifold with boundary. Generically, the partition function sums over branched covers. 
The derivative of the map at some points on the worldsheet can vanish, and the images of 
such points are called branch points. A particularly simple situation occurs when the Yang 
Mills theory is defined on a cylinder, and the area is taken to zero. Then the partition 
function sums only over unbranched covers. 

The large N expansion of 2dYM can also be expressed in terms of lattice gauge the- 
ory with S n gauge group which can be formulated following Wilson and letting edge 
holonomies take values in S n , with an appropriate plaquette action. This theory is topo- 
logical and will be called S n TFT (topological field theory). So along with the original 
U(N) description, there is the picture of S n TFT, as well as the Hurwitz space descrip- 
tion. The Hurwitz space connects most directly to the equations formulated on the string 
worldsheet, since the holomorphy condition is the localization locus for the string path 
integral. Explicit worldshseet actions have been proposed [6J [Tl El [9], [TO] . The 2d S n TFT 
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can be viewed as the string field theory of this Hurwitz space string theory. In this paper 
we will find the 2d S n TFT particularly useful. 

Other examples where the string-permutation connection is important is Matrix string 
theory[T7]. An orbifold superconformal field theory with target space S N (M. 8 ) arises as 
an IR limit of N = 8 SYM in two dimensions. Conjugacy classes in Sn correspond to 
string winding numbers which are dual to momenta of gravitons in the Matrix Theory 
interpretation. Yet another example is the connection between Belyi maps (a special class 
of branched covers with sphere target space) and correlators of the Gaussian Hermitian 
matrix model [TBI [T9] . 

From the above we take the lesson that the connection between permutations and 
strings is rather generic. In many cases, it is a tool to implement the ideas of strings 
emerging from the large N expansion in quantum field theories. Now a large number of 
counting problems can be formulated in terms of permutations acting on various types 
of sets. Some of these counting problems, such as the ones relevant to 2d Yang Mills 
and to Belyi maps, involve permutations acting on each other, for example by group 
multiplication. It will not surprise anybody that the question of counting Feynman graphs 
in quantum field theories, without large N, can be formulated in terms of symmetric groups 
acting on some sets having to do with vertices and edges. 

Investigating the counting of Feynman graphs more carefully reveals that, in fact, 
these counting problems can be formulated in terms of rather intrinsic properties of the 
symmetric groups themselves. The nature of the specific interactions determine the form 
of the vertices, and in turn certain subgroups of symmetric groups, which can be viewed as 
the symmetry of all the vertices. Then we draw in an idea from the study of Belyi maps 
and their associated ribbon graphs, which is called cleaning. This is the construction 
of associating to a given graph, a new graph obtained by introducing a new type of 
vertex in the middle of the existing edges, dividing them into half-edges. This allows the 
description of Wick contractions in terms of permutations which are pairings. This was 
exploited recently in a physics setting in [T9] . 

This line of thinking culminates in an elegant way to compute numbers of Feynman 
graphs of real scalar fields, which coincides with the results in a classic in the mathematical 
literature on graph counting [20]. Our derivation does not rely on the notion of graph 
superposition used in [20] and appears more direct, at least when approaching these 
graphs from the perspective of Wick contractions in perturbative quantum field theory. 
We also show that the formulae counting Feynman graphs, as well as those describing 
symmetry factors of individual graphs, can be interpreted as observables in the 2d S n 
TFT. The number n will depend on the number of vertices and edges in the Feynman 
graph calculation. We extend our considerations to QED, deriving new counting results 
on Feynman graphs. Again we obtain an expression of Feynman graph counting problems 
in terms of observables in the 2d S n TFT. An interesting aside is a surprising connection 
between the counting of QED Feynman graphs and the counting of ribbon graphs. 
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While the string-permutation connections underlie the emergence of strings from large 
N 2d YM, our results suggest that these same connections could also lead to the emergence 
of strings from quantum field theories such as that of a scalar field or QED. 

We now give an outline of the paper. In section 2, we review the connection between 
strings and permutations, which plays a prominent role in the string theory of large N 2d 
Yang Mills theory with U(N) gauge group. We describe the perspective of lattice gauge 
theory with S n gauge group and explain its topological nature by drawing on existing 
2dYM literature. We also review the Hermitian matrix model and its connection to Belyi 
maps, as another realization of the string-permutations connection. Finally, we mention 
the Burnside Lemma from combinatorics, which is used very generally for counting orbits 
of group actions. In section 3, we review Feynman rules for 4 theory. We recall how 
the symmetry factors get multiplied with additional group theory factors associated with 
global or gauge symmetry groups, and space-time integrals. Our main focus will be the 
combinatorics of the graphs and their symmetry factors. 

In section HJ we describe a method of enumerating Feynman graphs and calculating 
their symmetry factors, which relies on a pair of combinatoric objects. The first of the 
pair, which we call Eo, captures the form of the vertices. The second piece of data, which 
we call Ei, is associated with the Wick contractions. Given a graph, this data can be 
constructed by putting a new vertex in the middle of each edge. We can imagine these 
vertices to be coloured differently from the ones already present in the graph. This splits 
each existing edge into a pair of half-edges. We associate labels {1, 2, • • ■ } with each of 
these half-edges. This is explained in the context of real scalar field theory, in the first 
instance, with vacuum graphs in 4 theory. This formulation leads to simple one-liners in 
GAP (mathematical software of choice for group theory) [22] for calculations of symmetry 
factors and enumeration. This is explained with examples in Appendix O 

In section El we use the Burnside Lemma to count Feynman graphs, and obtain the 
first hints of stringy geometry of maps to a torus. The picture of maps to a torus has some 
intricacies, and it turns out that a deeper understanding of the combinatorics leads to a 
simpler picture in terms of strings covering a cylinder. This is is achieved by first recalling 
from the maths literature that graph counting formulae are expressed in terms of a certain 
double coset [20]. We find that the formulation of graph counting in terms of the pair 
(E , Ei) finds a natural meaning in terms of the double coset. This allows us to exhibit 
the counting of Feynman graphs in scalar field theory as an observable in 2d S n TFT on 
a cylinder. Following standard constructions in covering space theory, the S n data is used 
to construct unbranched covers of the cylinder. Some of the counting formulae that follow 
from the double coset picture can be understood using the idea of Fourier transformation 
on symmetric groups. Calculations explaining this are in Appendix [Bl Interestingly very 
similar calculations are relevant to recent results in correlators of BPS operators in A = 4 
super- Yang-Mills [2T] . 

We pursue, following [20], the application of the double coset picture to derive gener- 
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ating functions for graph counting in section |6j We also explain here how the formulae 
change when we generalize our considerations from 4 to 3 and other interactions. Gen- 
eralizations away from vacuum Feynman graphs to the case with external edges is also 
explained. The double coset formulae lead most directly to the set of all Feynman graphs, 
included disconnected ones. Counting formulae for the connected ones are obtained here. 
Appendix [D] contains explicit lists for Feynman graph counting. The first few terms agree 
with existing physics literature. A few of the series we consider are listed in the Online 
Encyclopaedia of Integer sequences [23], but the majority are not listed there. 

In section we show that these approaches work in a simple way for the ribbon 
graphs which arise in doing the 1/N expansion. While the ribbon graphs for large N are 
traditionally drawn using double lines, they can also be defined in a way closer to ordinary 
graphs, with single lines but with the difference that the vertices have a cylic order (see 
e.g. [21] [22])- This allows cyclic symmetry of the half-edges at the vertices in testing 
equivalence of differently labelled ribbon graphs, but not arbitrary permutations. Taking 
this into account, we express the question of how many Feynman diagrams correspond to 
the same ribbon graph (embedded graph) in group theoretic terms, associated with the 
symmetry breaking from permutation group to cyclic group. 

In section [HJ we apply these ideas to QED or Yukawa theory, giving the connection 
to 2d S n TFT and deriving generating functions. In section [HI we adapt the counting to 
QED, with the vanishings due to Furry's theorem taken into account. Again we obtain 
the 2d S n TFT connection and the generating functions. Manipulations of the double 
coset description for QED leads to a somewhat surprising connection between QED graph 
counting and ribbon graph counting. We explain this by describing a bijection between 
these graphs. The arrows on electron loops provide the cyclic symmetry which is key to 
ribbon graphs. A consequence is that the number of QED/Yukawa vacuum graphs with 
2v vertices is equal to the number of ribbon graphs with v edges. 

In section \TU\ we discuss our results with a view to extracting some lessons for gauge 
string duality. Ribbon graphs give a vivid physical picture of how strings arise from 
quantum field theory, but an equally compelling and arguably simpler physical picture is 
that of strings winding on a circle being described by permutations. This latter picture 
has been exploited here to give new connections between QFT Feynman diagrams and 
string counting. The natural question is whether this connection can be extended to 
extract from Feynman graphs of QFT, without large N, something more than stringy 
combinatorics to include space-time dependences of correlators and S-matrices. 

In section [TTJ we summarize the paper and discuss some concrete avenues for future 
work. Appendix [A] describes a semi-direct product structure of the Automorphism groups, 
whose orders give symmetry factors of Feynman graphs, and points out a difference be- 
tween the notion of Automorphism most commonly used in the mathematics literature 
on graphs and the one relevant to symnmetry factors of Feynman graphs. 
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2 Review of strings and permutations 



Consider a string wrapping a circle. The winding number is a topological characteristic 
of the map. In string theory, the string is part of the worldsheet, the circle is viewed as 
part of target space. Let the target circle be parametrized by X with X ~ X + 2n. Let 
the string be parametrized by a with a ~ a + 2-k. A string with winding number k is 
described by 



X = ka 



(2.1) 



For multiple strings, we may have distinct winding numbers. Adding the winding numbers 
gives the degree of the map from the strings. 

Given such a configuration of winding strings, we may label the inverse images of a 
fixed point on the circle from 1 to n. Following the inverse images as we move round the 
target circle yields a permutation. The two possible winding configurations at n = 2 are 
shown in Figure HJ 




(1)(2) 




(12) 



Figure 1: Winding strings and permutations 



In ordinary string theory, the winding number is one of the quantum numbers speci- 
fying a string state. In some topological settings, the string winding number is the only 
relevant quantum number. 



2.1 Review of results from 2d Yang Mills 

Recall that 2d Yang Mills on a target space with a boundary has the partition function 

Z = J2(VimR) 2 - 2G - B e- Ac ^ Yl Xr{Ui) (2.2) 

R i 

G is the genus of the surface, B is the number of boundaries, f/j are the holonomies at 
the boundaries and A is the area of the target and C*2(-R) is the quadratic Casimir [1]. 
The choice of fixed boundary holonomies is illustrated in Figure [2j 
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u. 



Figure 2: The 2dYM partition function with boundary holonomies is known exactly 



For a target cylinder, we have 

Z(U lt Ui) = 5>*(£/i)x«(£/ 2 )e- AC2 ^ 



(2.3) 



R 



in terms of characters of U(N) elements. The power of the dimension is zero because 
G — 0, B — 2. The characters have an expansion in multi-traces. These traces can be 
constructed from permutations as follows 



tr(aU) =UV- W 1 ■■■U i i 

y ' * CT (1) V(2) 1 



II 

tr(n) 



(2.4) 



This is a very useful formula when developing a string interpretation for Wilson loops in 
2dYM [26]. 



Using these observables for the cylinder, we define 



Z(<7i,<7 2 ) = j dUKlUitritnUbtr&VtWUuUz 



This integral can be done, and in the zero area limit A = we find 



(2.5) 



(2.6) 



Rhn 



The notation R\- n means that R is running over partitions of n (row lengths of Young 
diagrams) which parametrize irreducible representations of S n , and the X-r( (J ) are char- 
acters of S n elements. Using the invariance of the character under conjugation, we can 
replace 



Xr< 



0-2) = ^7 V ^(7^27 l ) 



(2.7) 
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This gives 



Z(o-i,0- 2 ) = — ^ ^Xr{°i)Xr(ip21 X ) (2.8) 

' jeSn Rrn 

Since ^ 7 7027 -1 is a central element in the group algebra of S n , we can use Schur's 
Lemma to combine the product of characters into a single character 

Z(ai, cr 2 ) = — ^flXflC^iT^T -1 ) ( 2 - 9 ) 

' 7G5n R\-n 

Now use the Fourier expansion on the symmetric group, which allows the delta function 
to be written in terms of characters, to obtain 

£>\Pu °~2) 

7 

When the gauge group is a product U(N) x U(N) x U(N), we have to specify a 
holonomy for each gauge group in the product. Denote them by U, V, W. Recalling that 
the character of a direct product is the product of the characters we have the zero area 
partition function for a cylinder 

Z{Ui, Vt, Wi; U 2 , V 2 , W 2 ) = Xr(Ui)xs(V 1 )xt(W 1 )xr(U 2 )xs(V 2 )xt(W 2 )) (2.11) 

R,S,T 

The boundary observables in this case, in a basis appropriate for a string interpretation, 
are tr(aU)tr(pV)tr(rW). In terms of permutations 

Z(a 1 ,p 1 ,r 1 ;a 2 ,p 2 ,r 2 ) = ^ ° ft ° r i 7 ^ ° ft ° ^2 T' 1 ) (2-12) 

We have discussed above the classification of string maps to a circle target space, in 
motivating the choice of boundary observables in large N 2dYM. The full interpretation 
of this result requires an extension of our considerations to a Riemann surface target. In 
this case, the strings-permutations connection has deeper implications captured in the 
Riemann existence theorem. Consider the equivalence classes of holomorphic maps from 
the worldsheet (E/J to the target (£(?) with two maps fi and f 2 defined to be equiva- 
lent if these exists a holomorphic invertible map <p : — > such that fi = f 2 ° 
Given a holomorphic map (branched cover) with L branch points and of degree n, we 
can obtain a combinatoric description by picking a generic base point and labeling the 
inverse images as {1, 2, • ■ ■ ,n}. By following the inverse images of a closed path starting 
at the base point and encircling each branch point, we can get a sequence o~\, o~ 2 , ...,o~l of 
permutations. For G > 0, there are also permutations for the a, b cycles of Eg, denoted 
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Si,ti for i = 1 ■ • ■ G. Two equivalent holomorphic maps are described by permutations 
<7i, cr 2 , (Jl, Si, ti, • • • , s<3, to and c^, cr^,, cr^, • • ■ , s' G , t' G which are related by con- 

jugation (jj = acr-a -1 ,^ = as[a~ l ,ti = at[a~ l . This correspondence between sequences 
of permutations and holomorphic maps is captured in the Riemann existence theorem (see 
for example [2S])- Relations in the fundamental group of the Riemann surface punctured 
at the branch points are reflected in the sequence of permutations. 

The delta function in (12.101) enforces a relation in the fundamental group of the cylin- 
der. When this partition function is generalized beyond zero area, the sum fl 2 . 1 [) is 
modified to include additional permutations which can be interpreted as a counting of 
branched covers. Even at zero area, but for other Riemann surfaces, there are additional 
permutations associated with branch points. We will come back to this point in the 
discussion of (I2.16p . 

2.2 Lattice topological field theory for S n 

In lattice gauge theory [2T] for two dimensions, we discretize (e.g triangulate or give a 
more general cell decomposition for) the Riemann surface. To each edge we associate a 
group element. To each plaquette, we associate a weight which depends on the product of 
group elements along the boundary of the plaquette. Let us choose the plaquette weight 
to be 



where a is the product of group elements around the boundary of the plaquette. The 
partition function for the discretized manifold is the product of weights for all the plaque- 
ttes. The partition function can be shown to be invariant under refinement of the lattice 
so that the result depends only on the topology of the Riemann surface. This is proved 
in [28]. The language used there is that of continuous groups, but by replacing integrals 
with sums, it applies equally well to finite groups. In [28] surfaces with a finite area A 
are considered. For the topological aspects of the theory that interest us, it is enough 
to consider A = 0. For manifolds with boundary one fixes the group elements at the 
boundaries. As a simple example, consider a lattice theory defined on a disc. The target 
might be discretized with either one plaquette or two plaquettes as shown in Figure [3j 
Starting from the discretization given in (a) we find 




(2.13) 



R\-n 



Z(a 1 ,a 2 ) = ^2,Zp{ai^{ 1 )Z P (aj) 
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7 



5{a l a 2 ) 



(2.14) 
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which is precisely the partition function obtained using the discretization shown in (b). 
The generalization to finer discretization works in exactly the same way. 



(a) 





Figure 3: Two discretizations of the disc with the same boundary condition. 



For the cylinder 

= X)5(<W27 _1 ) (2-15) 



which is of course (I2.10p . For a 3-holed sphere 

^720-372 



[n\] 

7l>72 



-T 



Hl R dl 



(2.16) 



This is also derived in [29J from a more axiomatic approach. 

These formulae from lattice gauge theory of S n arise from the leading the large N limit 
of the U (N) gauge theory at zero area. The cylinder case is special in that this is the full 
answer to all orders in the 1/N expansion. The answer in (12 . 1 6[) is obtained after using 
the leading large approximation of Q, which is a sum over elements in S n weighted by 
powers of N. See for example 6.1.2 of [6] for a discussion of Q. At subleading orders, the 
Q factor contains a sum over permutations, which can be interpreted in terms of sources 
of curvature in the S n bundle. From the Hurwitz space string interpretation of 2dYM, 
there are branch points in the interior. 

The S n lattice gauge theory perspective for U (N) gauge theory at large iV is empha- 
sized in [30J. Computations of 2dYM partition functions for general Riemann surfaces 
with boundaries expressed in the symmetric group basis, and the connection with Hurwitz 
space counting is given in [31 EJ [26] . The connection to topological field theory with S n 
gauge group was also observed in [3T] . 
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To summarize, the large N expansion of 2dYM with U(N) gauge group on a Riemann 
surface £g,b (genus G, with B boundaries) can be expressed in terms of symmetric groups. 
This combinatoric data arising has an interpretation as 2d gauge theory with S n gauge 
group. This connects directly with the string theory interpretation in terms of maps from 
a string worldsheet (genus h related to the order in the 1/N expansion) using classic 
results relating the space of branched covers (Hurwitz space) to symmetric groups. 

We will find in this paper that counting problems of Feynman graphs can be expressed 
in terms of certain generalizations of the 2dYM results, which can be expressed elegantly 
in terms of the 2d S n gauge theory perspective, and also admit a connection to branched 
covers. 

2.3 Permutations and Strings beyond 2dYM 

The connection between strings and permutations also has interesting implications for the 
correlators of the Gaussian hermitian one-Matrix model [TH]. The computation of corre- 
lators can be mapped to the counting of certain triples of permutations which multiply to 
the identity. These count holomorphic maps from world-sheet to sphere target with three 
branch points on the target. Holomorphic maps with three branch points are related, 
by Belyis theorem, to curves and maps defined over algebraic numbers. Thus, the string 
theory dual of the one-matrix model at generic couplings has worldsheets defined over the 
algebraic numbers and a sphere target. For related ideas see [32j [33] . 

Finally, yet another connection between strings and permutations that is relevant to 
our present discussion is provided by Matrix String Theory, defined by the IR limit of 
two-dimensional N = 8 SYM. This limit is strongly coupled and a nontrivial confor- 
mal field theory describes the IR fixed point. The conformal field theory is the N = 8 
supersymmetric sigma model on the orbifold target space (R 8 ) N / Sn, formed from the 
eigenvalues of the Higgs fields Xj of the theory [T7]. If we go around the space- like S 1 of 
the world-sheet, the eigenvalues can be interchanged. Concretely 

X I (a + 2ir)=X g{I) (a) 

where g takes value in the Weyl group of U(N) which is the symmetric group Sn- These 
twisted sectors correspond to configurations of strings with various lengths. Consequently, 
twisted sectors with given winding number are labeled by the conjugacy classes of the 
orbifold group Sn- 

2.4 A useful theorem in combinatorics 

Let G be a finite group that acts on a set X. For each g in G let X 9 denote the set of 
elements in X that are fixed by g. Burnside's lemma asserts the following formula for the 
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number of orbits 

Number of orbits of the G-action on X = — ^ \ xg \- ( 2 - 17 ) 

I I g&G 

Thus the number of orbits is equal to the average number of points fixed by an element 
of G. This is called the Burnside Lemma or the Burnside counting theorem. Useful 
references for the Burnside Lemma are [JT] [12], the former also provides other useful 
combinatoric background. 



3 Review of Feynman graph combinatorics 

There are many excellent articles and textbooks that deal with the Feynman graph ex- 
pansion of perturbative quantum field theory. Here we will simply review those aspects 
most relevant for our goals. For a relevant reference that has more details see |43j . 

Perturbative quantum field theory expresses quantities of interest (for example, an 
amplitude A) as a power series expansion in the small coupling g 

oo 

A(g) = J2^9 v (3.1) 

fc=0 

The coefficients A v are obtained by summing Feynman graphs with v vertices D v 

A v = J2 C d,N Dv F Dv (3.2) 

D v 

Cd v is the symmetry factor of graph D v , Nd v is a group theory factor coming from 
the color combinatorics of global or gauge symmetry groups and Fjj v is the result of 
integrating over the loop momenta in D v . For concreteness assume that g is the strength 
of a 4 interaction. In this case the factor Np v is 1. If we canonically quantize the 
theory, we expand about the free theory using Wick's Theorem. The Feynman graphs 
are used to compute the sum over all possible Wick contractions. Different ways of doing 
the Wick contractions can lead to the same Feynman graph. The symmetry factors Cn v 
keep track of this. Perturbation theory expands the exponential of the g§ A interaction. 
Consequently, g v come with a ljv\ from the Taylor series of the exponential. Accounting 
for this factor, the number of Wick contractions leading to Feynman graph D v is v\Cd v - 
This observation can be used to generate an interesting sum rule for the symmetry factors. 
The sum of graphs with v vertices and E = 2n external lines reproduces the complete set 
of Wick contractions for Av + 2n fields. The total number of Wick contractions is equal 
to the number of pairs that can be formed from the Av + 2n fields and is also equal to the 
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sum over D v of v\Cr> v . Consequently 

v] 52 c d v = (4« + 2n-l)!! (3.3) 

D v 

Note that for this sum rule to hold we must not subtract vacuum graphs. This sum rule 
provides a good test to check the symmetry factors. One way to count how many Wick 
contractions will produce the same graph, is to count the number of ways of interchanging 
components which don't modify the graph. These are the automorphisms of the graph. 
Indeed, the standard recipe jS] for computing the Cn k starts by giving a set of rules to 
determine these automorphisms. For g0 4 theory the rules are 



• for a closed propagator (one with both ends attached to the same vertex) there is a 
swap which exchanges the endpoints. 

• for p propagators which each start at the same vertex and end at the same vertex 
(the starting and ending points might not be distinct) there are p\ transformations 
that permute the propagators. 

• for n vertices that have the same structure there are n\ transformations that permute 
the vertices. 

• for d identical graphs there are d\ transformations that permute the graphs. 



Denote the automorphism group of a Feynman graph D by Aut(-D). The order 
|Aut(£))| is the product of the numbers obtained applying each rule above. A Feynman 
graph built using v vertices comes with a coefficient 

c °- = j^m < 3 - 4 > 



Applying the rules to compute symmetry factors can be tricky. The excellent textbook [ 
suggests that "When in doubt, you can always determine the symmetry factor by counting 
equivalent contractions" . As v grows, this quickly becomes hopeless. One of the results of 
this paper is give conceptually simple permutation group algorithms for symmetry factors 
(Appendix [C]). 

One basic quantity of interest is the number of graphs 

N V = J2 X ( 3 - 5 ) 

D v 

contributing, since it is an important effect in determining the behavior of field theory 
at large orders in perturbation theory [45] . For questions of this type the value of Fjj v is 
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unimportant (we are implicitly assuming they are all roughly the same size) so that we 
may, for simplicity, work in zero dimensions which amounts to setting Fjj v = 1. 

Instead of focusing the discussion on any particular amplitude, it is useful and stan- 
dard, to study the generating functional of correlation functions 

Z[J]= J [D<p\e lS+i $ dixJ4, (3.6) 

where 

S = j d 4 x Qfy ^ " \™ 2 <P 2 ~ W> 4 ) (3-7) 

The value of Z at J = computes the sum of all vacuum graphs. Derivatives of Z with 
respect to J compute the sum of all Feynman graphs with E external legs 



1 () n 7 \ 71 

(Q\ (j) ( Xl ) ( f ) ( X2 )... (j) ( XE )\0} = - 



i E SJ(xi)5J(x 2 ) ■ ■ -5J(xe) 



(3.8) 

j=o 



Another quantity of interest is the logarithm of Z, W = InZ. The value of W at J = 
computes the sum of all connected vacuum graphs. Derivatives of W with respect to 
J, at J = 0, compute the sum of all connected Feynman graphs with E external legs. 
This suggests another interesting question: how many connected Feynman graphs iV™ nn 
are there? The connection between Z and W holds for graphs with the symmetry factor 
Cd v included. The relation between iV™ nn and N v (which does not include the symmetry 
factors) is different. 



4 Feynman graphs in terms of the pair (Eo, E] 

Consider a vacuum Feynman graph in <p A theory. Let v be the number of vertices. At 
v = 1 we have a single graph, given in figure HI At v = 2 we have three graphs given 



Figure 4: One vertex vacuum graph in 4 theory 



in figure [5j To get a systematic counting, we describe these graphs in terms of some 
numbers. One way to do this is to label the vertices 1, 2, • • • ,v, and list the edges as [ij]. 
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ooo 




Figure 5: Two vertex vacuum graphs in 4 theory 



This method is well-known in graph theory and quickly leads to elegant results for graphs 
where there are no edges connecting a vertex to itself, and no multiple edges between a 
given pair of vertices. In the case at hand, we do have loops and we do have multiple 
edges between the same vertices. For this reason, the graphs at hand are sometimes called 
multi-graphs in the mathematics literature. 

For the graphs we consider, it is more convenient to give a combinatoric description, 
by putting a new vertex in the middle of each edge. Each edge is thus divided into two 
half-edges. Each half-edge is labeled by a number from 1,2, • • -Av. Using this labeling, 
the graphs shown in figure H] and figure [5] have been labeled in figure [6] and figure [7J 
respectively. We have drawn the original vertices as black dots and the new vertices as 
white dots. 




Figure 6: Numbering the half-edges 



Now the graphs can be described by specifying the list of 4-tuples of numbers at the 
black vertices and the pairs of numbers at the white vertices. We thus introduce two 
quantities E , Si. For the v — 1 graph, with the labeling chosen in figure [6] 

5 = <1,2,3,4> 

51 = (12)(34) (4.1) 

For the three graphs at v — 2, inspection of figure [7] shows that we have fixed 

(a) S = < 1,2,3,4 >< 5,6,7,8 > 

Si = (12)(34)(56)(78) 

(6) S = < 1,2,3,4 >< 5,6,7,8 > 
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(a) 



1(J2 5(J6 
3(j4 71)8 



2 1^6 5 

(b)OOQ 

3 4 7 8 




Figure 7: Numbering the half-edges 



(c) 



Ei = (23) (16) (47) (58) 

S = < 1,2,3,4 >< 5,6,7,^ 
Ei = (15)(26)(37)(48) 



> 



(4.2) 



Clearly, we could have labeled things differently producing different pairs Eo, Ei. These 
relabelings are elements of 5*4 in the case v — 1 and elements of Ss in the case v — 2. 
More generally, we have Si v . It is useful to consider breaking the half-edges, leaving us 
with 4t> half-edges dangling from v 4-valent vertices. And 4v half-edges from 2v 2-valent 
white vertices. We imagine we did not remove the labels as we broke the half-edges, so 
we know how to put the graph back together: by connecting half-edges with identical 
numbers. After breaking the half edges in figures [6] and [7] we obtain figures [8] and [9l 



2 
4 



y\ A 

1 2 3 4 



Figure 8: Splitting the half-edges 



The symmetry group of the 4-valent vertices includes the 4! permutations of each 
vertex, along with the v\ permutations of the vertices themselves. These permutations 
form a subgroup of S± v called the wreath product S V [S^\. It is the semi-direct product 
of (S4)" x S v . The physics reader is not required to have any prior knowledge of these 
groups. All the relevant facts we will need from the math literature will be quoted as we 
need them. 

The symmetry group of the Feynman graph, whose order appears in the denominator 
as the symmetry factor, is obtained from the action of the symmetric group S^ v on the 
pairs (Eo, Si). A permutation in S^ v is a map from integers {1, 2, • • -4w} to these integers 
{1,2, • • ■ , Av}, the map being one-one and onto. It can be written as an ordered list of 
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( a ) X X 

( b ) X X 

Figure 9: Splitting the half-edges 

the images of {1, 2, • • • , 4v}, e.g a cyclic permutation of {1, 2, 3, 4} can be written as 

2341 (4.3) 

Alternatively we use cycle notation, (1234) for the example at hand. The image of every 
integer in the bracket is the one to the right, the image of the last is the first. 

A permutation a acts on E , Si by taking the integers i in these tuples to cr{i). We 
denote this action by 

^(Eo.EO^HEo),^!)) (4.4) 

If (E' , Tii) = (cr(So), <t(£i)) for some a, then both pairs (E^E^) and (Eo, Si) represent 
the same Feynman graph. We may say that Feynman graphs correspond to orbits 
of the S4,, action on the pairs (E ,Ei). 

As we saw in the examples above we can fix the form of Eo- The subgroup of S± v 
which preserves this is the S^Sy. Then different Feynman graphs correspond to orbits 
of S^Sy acting on the pairings Ei. These are nothing but elements of the conjugacy 
class [2 2v ] in S^ v . This leads directly to the computation of numbers of Feynman graphs 
and symmetry factors using GAP software, which is the mathematical package for group 
theory. See Appendix O 

5 (j) A theory and string theory 

In this section we develop a number of key ideas. As we saw in the last section, to generate 
all possible Feynman graphs, we can fix Eo and allow Ei to vary. Using the action 

c : (Eo, Ex) ( C 7(E ),a(E 1 )) (5.1) 

of er g 5*4„, the condition cr(E ) = E fixes a to be in the subgroup 5 l ,[5 , 4 ] of S± v . Per- 
mutations Ei that can not be related by a G S^SV] are distinct Feynman graphs. This 



A A A A 

1 2 3 4 5 6 7 8 

A A A A 

2 3 1 6 4 7 5 8 

A A A A 

1 5 2 6 3 7 4 8 
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leads to the realization of Feynman graphs as orbits of the vertex symmetry group acting 
on Wick contractions. A simple application of the Burnside Lemma in section 15.11 then 
gives an explicit formula (15. 4 p for the number of Feynman graphs. The formula has an 
immediate interpretation as a sum over covers of a torus. This suggests an interpretation 
as string worldsheets mapping to a torus. However, the sum over maps to torus have some 
constraints which seem non-trivial to implement geometrically. This motivates a deeper 
look at the combinatorics, which reveals a somewhat simpler connection to strings with 
a cylinder target space. In section 15.21 we explain that Feynman graphs are elements of 
a double coset. This description is already present in [2U], where it is derived using an 
operation called graph superposition. Our derivation does not use this operation and is 
obtained more directly using the half-edges introduced in Section HJ The pair (Eo, Si) 
encountered there finds a natural interpretation in terms of the double coset. This is a 
key result with a number of implications. It allows us to express the counting of Feynman 
graphs as well as their symmetry factors in terms of data in S n TFT. In turn this S n 
TFT data is used to construct covers of the cylinder. The figure [TU] summarizes the key 
message of this section. This will set the stage for Section where generating functions 
and counting sequences will be given. 



(|) Theory 



Vacuum graphs with 
v vertices 



H A S 4vM 



Double Coset 



String Theory 



ri 



Figure 10: Double coset connection 



5.1 Commuting permutations 

A vacuum Feynman graph in 4 theory is constructed by starting with v vertices, each 
having 4 edges attached, and connecting up the edges, as shown in Figure 11. Without 
any loss of generality label the edges of the first vertex (1, 2, 3, 4), the edges of the second 
vertex as (5, 6, 7, 8), etc. We could have chosen different numbers (from (l.Av) to go with 
each vertex), but we have made a choice. This is a choice of 

E = (1,2,3,4>(5,6,7,8> ■■■ (5.2) 
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where the angled brackets are completely symmetric and the different brackets can be 
freely interchanged. The symmetry of E is S v [S4] 

To complete the diagram we need to give the Wick contractions. The Wick contrac- 
tions are specified by choosing an element Ei in [2 2v ]. For example 

(12)(34)..(4u-l,4v) (5.3) 

corresponds to the disconnected graph with a bunch of eights. 

To construct all possible Feynman graphs we need to allow Ei to run over all possible 
Wick contractions holding E fixed to the choice we made above. This is easily done 
by acting on Ei with S^Sy. It follows that the distinct Feynman graphs correspond to 
orbits of the vertex symmetry group S^fSy acting in the set of pairings [2 2v ], which are 
elements of a conjugacy class in S^ v . If we apply the Burnside Lemma, we see that the 
number of Feynman graphs is equal to 

Number of Feynman Graphs = - 5(7<T7 _1 cr _1 ) (5-4) 

This is a sum over certain covers of the torus of degree 4v. The worldsheet is also a torus. 
One of the monodromies is constrained to lie in the subgroup Hi = S V [S^], the other lies 
in the conjugacy class C = [2 2v ]. Constrained sums of this sort where the permutations 
are constrained to lie in a given conjugacy class are naturally motivated from holomorphic 
maps. Here we are constraining one permutation to lie in a subgroup and one to lie in a 
conjugacy class. Such a constraint is entirely possible in S n TFT, but the the geometrical 
interpretation in terms of Hurwitz space is, at this point, less clear. Subsequent formulae 
we derive will show the emergence of a cylinder rather than a torus. The formulae then 
become more symmetric with respect to exchange of H\ and Hi- 

For purposes of counting applications, the equation (15 .4p is easy to implement in GAP, 
although related techniques involving cycle indices introduced in Section [6] will be more 
efficient. 

Apart from the number of Feynman graphs, it is also interesting to compute the 
automorphism group of a given graph. This appears, for example, in the denominator 
for the formula H3.4[) for the symmetry factor. With our choice of a fixed Eo, a graph is 
uniquely specified by E^ For this reason we will denote the automorphism group of a 
given graph by Aut([E 1 ]// 1 ). The [E^^ denotes an equivalence class under conjugation by 
H\. The automorphism group is given by those elements of Hi which leave E x invariant. 
The order of the automorphism group is 

lAutCpiU)! = ^il' 1 ^ 1 ) (5.5) 

76^1 

1 Some authors prefer the notation S4 I S v . We follow [20! for the reason he gives : that the S^Sy 
notation connects nicely with the substitution formula for cycle indices 
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A symmetry 7 and a Wick contraction Si specify an unbranched cover of the torus, of 
degree 4t> with 7 constrained to lie in Hi and Si constrained to be made of 2-cycles. 

Recall from section [3] that the number of Wick contractions leading to a particular 
Feynman graph is 

= (5 ' 6) 

Noting that Aut(S ) = (S 4 ) v x S v we find 

\(S,Y\ [Aut(E )l 
'lAutCpilffJI |Aut(£>)| 1 ' J 

This formula has a very natural interpretation: start by making a choice of S . To deter- 
mine a Wick contraction, we need to specify the cycle £1. Any other Wick contraction 
contributing to the same graph can be obtained by swapping vertices or swapping the half 
edges at a given vertex. The full set of these transformations is performed by Aut(So). 
If we have an element of Aut(So) that leaves Si invariant (this by definition belongs to 
Aut(-D)) we do not get a new Wick contraction. Thus, the formula (15.61) is an application 
of the Orbit-stablizer theorem 



5.2 Double Cosets and the meaning of the pair (So, Si). 

Read [20] derives explicit counting formulae for graphs by developing a formulation in 
terms of double cosets. He arrives at the double cosets by a procedure of graph super- 
position. We will arrive at the same double coset descriptions using the procedure of 
separating the edges into half-edges, and keeping track of the permutations which contain 
the information of how the half-edges are glued to make up the graph (See Figure [9]). 

There are two related double coset descriptions relevant to 4 graphs. The first one is 
a double coset of S n x S n . The second is a double coset of S n and can be obtained by a 
gauge fixing of the S n x S n picture. Let us start with this description. The elements of 
the double coset are pairs 

(cri,(7 2 )e(S n xS n ) (5.8) 

with the equivalence 

(a u a 2 ) ~ (aa 1 (3 1 ,aa 2 (3 2 ) (5.9) 

where n = 4t>, a G S iv , f}\ 6 H\ — S v [Sa\ and j3 2 e H 2 = S 2v [S 2 }. 

To understand why this double coset counts Feynman graphs of </> 4 theory consider 
graphs with v 4-valent vertices which will have a So, Si description 

(S ,Si)G(<4^>,[2 2 ^) (5.10) 
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As explained in section H] the graph splits up into 2v bi-valent white vertices and v 4- 
valent black vertices. Label the edges connected to the white vertices 1 ■ ■ • 4v, so we have 
edges {Wi, W 2 ■ ■ ■ W iv }. Label the edges connected to the black vertices {1, 2, • • ■ , 4t>}, 
so we have {Bi, B 2 , ■ ■ ■ B^ v }. All possible relabellings of the edges of the white vertices 
are paramatrized by o\ G S± v . All possible relabellings of the edges of the 4-valent black 
vertices are parametrized by a 2 G S^. Given any (01,02) we can construct a graph by 
gluing 

W Ul(i) B a2{i) (5.11) 

for i — 1 • ■ • 4v . 

Clearly by considering all possible (01,02) we can get all possible graphs. 

If we replace the % by a(i), we get the same graph, since the labelings do not matter. So 
we learn that (aa±, aa 2 ) and (aai, ca 2 ) produce the same graph. Hence we are interested 
in equivalence classes 

(<7i,<r 2 ) ~ (<wi, 0:0-2) (5.12) 



We also know that the disconnected graph of bi-valent vertices has a symmetry of 
H 2 = S 2v [S 2 ], so that cr 2 and a 2 f3 2 with f3 2 G H 2 is the same graph. In the same way, 
<7i ~ cri/3i with Pi G Hi. This leads us to the conclusion that Feynman graphs are 
equivalence classes of the equivalence relation 



(0-1,0-2) ~ (aoi/3i, aa 2 /3 2 



(5.13) 



completing the demonstration that Feynman graphs are in one-to-one correspondence 
with elements of a double coset. 
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Figure 11: Double coset connection 



If we choose a = (3 t a x , then any pair is mapped to 

(l^r^rWa) (5-i4) 
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So now only the combination a 1 a 2 appears which is in S n , and we have equivalences by 
right multiplication with H 2 and left multiplication with Hi. We learn that 

S n \ (S n x S n )/{Hi x H 2 ) = H x \ S n /H 2 (5.15) 

From (I5.13P we can consider first doing the coset by Hi and H 2 , to we have an action of 
S n on S n /Hi x S n / H 2 . This is the description we are using when we impose S n relabeling 
equivalences on (£ , £1) which are nothing but a parametrization of S n /Hi x S n /H 2 . So 
we have now learnt the group theoretic meaning of the Eq,£i, which we wrote earlier 
as a combinatoric description of the graph that followed (at the beginning of section 4) 
immediately from introducing the white vertices to separate the existing vertices and keep 
track of the Wick contractions. 

This line of argument makes the generalization to ribbon graphs clear. Not any per- 
mutation of the legs connecting to a black vertex is allowed: to preserve the genus of the 
ribbon the cyclic order of the labeling must be respected. Consequently, Hi = S V [S^ 
must be replaced by Hi = S^Z^]. We pursue this direction in section [71 



5.3 Cycle Index formulae related to double cosets 

In this section we would like to make a connection with the classic results of Read [20] on 
the counting of locally restricted graphs. The starting point of [20] it to count equivalence 
classes of pairs (01,02) in (S n x S n ), where the pair (01,02) is equivalent to (61,62) if 

(61,62) = (xaigi,xa 2 g 2 ) (5.16) 

where x G S n and gi E Hi, g 2 e H 2 . This is equivalent to counting elements of the double 

coset 

S n \ (S n x S n )/(Hi x H 2 ) = Hi \ S n /H 2 (5.17) 

The number of equivalence classes is counted by N(Z(H{) * Z(H 2 )), where the Z's are the 
cycle indices, and N is a star-product for multivariable polynomials defined as follows: 
Consider two polynomials P and Q in variables j\, f 2 ..f n 

p — V* P ■ ■ P 1 P 2 ■ ■ ■ n 

jhn 

Q = zZQn,n-,Uif 3 2---Pn (5-18) 

jhn 

where j is a partition of n : ji + 2j 2 + • • • + nj n = n. We abbreviate the coefficients as 
Pj, Qj. The function N of a star product is now defined by 

N(P*Q) = Y / P j Q j Sym(j) (5.19) 
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where Sym(j) is the symmetry of the conjugacy class corresponding to partition j 

n 

Sym(j) = Hi/ 



i=l 



For our applications we will sometimes want to change the numerator group, in which case 
the symmetry of the conjugacy class will change. This generalization of Read's formula 
is discussed in section 18.11 

Read [20] also proves 





1 


\Hi 





NiZiHi) * Z(H 2 )) = Z^ ] Z^Sym(p) = 1 E u(a u a 2 ) (5.20) 

where 1/(01,02) is the order of the intersection aiHia^ 1 U a 2 if 2 G t 2~ 1 - We can rewrite the 
above as 

NiZm * Z(H 2 )) = ^rj^j E E E ^aiWW^ 1 ) 

Ol, OaES^v UlS-ffl U2&H2 

E E E ^kk 1 ^ 1 ) (5.21) 

feS54„ U1GH1 U2&H2 

which is reminiscent of a cylinder partition function for gauge theory with gauge 
symmetry, and with holonomies at the boundaries restricted being summed in Hi and if 2 
respectively. 

Since these formulae will come up, in slight variations, repeatedly, we make some 
comments on notation. We will denote the number of points in H\ \ G/H 2 as J\f(Hi, H 2 ) 
when G is just the symmetric group S n , or Af(Hi, H 2 ] G) more generally. The function of 
N of star products will also sometimes written to make the numerator group explicit. So 
the previous formulae can be expressed as 

M{H X , H 2 ) = N{Z{H X ) * Z{H 2 )) (5.22) 

or 

M{H U H 2 ; G = S n ) = N{Z{H l ) * Z(H 2 ); G = S n ) (5.23) 
Note that we can also write (I5.2ip as 

WT E E E Xr^Xr^) (5.24) 



\H\ \\H2\ 



RhSiv uiGHi U2GH2 



The sums over u\ and u 2 produce projection operators which project onto the trivial 
representation so that the only representations R which contribute are the representations 
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of S^v which contain the trivial of S^Sy and the trivial of S24S2]. Each such R contributes 
the product of the multiplicities with which the trivial of Hi and H 2 appear. 

Number of Feynman Graphs = ^ Xf fl -M-i H (5.25) 

R\-4v 

We have used the notation Mf H for the multiplicity of the one- dimensional representation 
of Hi when the irreducible representation R of is decomposed into representations of 
the subgroup Hi. 



5.4 Action of Vertex symmetry group on Wick contractions 

We now have two different ways to compute the number of Feynman graphs: as the 
number of orbits of the vertex symmetry group acting on Wick contractions, which leads 
to 

E E ^i^rv- 1 ) (5.26) 

or in terms of the star product of the cycle indices N (Z (Hi)* Z (H2)) . We will demonstrate, 
using some general group theory, the equivalence of (15. 26ft and (15.2ip . 

Noting that Hi = S V [S4], (15.261) and (15.211) are equivalent if 

Simau^a- 1 ) = J-r E KuxWb- 1 ) (5.27) 

a£[2 2v ] ' 2 ' b€S 4v u 2 €H 2 

To see that this relation holds, suppose there is a solution to the first delta function, a is 
conjugate by some b to ctq = (12)(34)..(4t> — l,4v) 

6~V 6 = a (5.28) 

We know that all elements commuting with <To are in if 2 = Si v [S2] ■ From ( 15.27ft we have 

uiou^ 1 = a (5.29) 

So 

uib^aobui 1 = rV 6 (5.30) 

It follows that 

bmb' 1 = u 2 (5.31) 
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for some u 2 in H 2 . The b that takes ctq to a can be multiplied by an arbitrary element of 
H 2 . So if we replace the sum over a G [2 2v ] by a sum over b, u 2 , we will be over counting 
by \H 2 \. So we conclude that the equality (15.271) is correct. 

It is instructive to note that the set of elements [2 2v ] can be thought as the set of 
cosets S4 V /H 2 , where H 2 = S 2v [S 2 }. To see this, start by noting that S iv acts on cr by 
conjugation to generate all the elements in [2 2v \. For any a G [2 2v ], we have 

a = a<Toa~ l (5.32) 

for some a G S^ v . If we multiply a on the right by j3 G S2 V [S 2 ], we get the same a. We 
can write 

[2 2v ] = S Av /H 2 (5.33) 

Thus, we can rewrite (I5.26P as 

Number of Feynman Graphs = . ^(uxcru^V -1 ) (5.34) 

^ '> V ' meffi aes iv /H 2 

A similar argument implies that 

< 4" >= Sto/H! (5.35) 

Here < 4 U > stands for the space of E 's. 

Using that fact that the expression (I5.2ip is symmetric under the exchange of Hi and 
if 2; we can also write 

Number of Feynman Graphs = ^ 5(tt 2 (E ), E ) (5.36) 

Thus, by Burnside's Lemma the number of Feynman graphs is also equal to the number of 
orbits of S 2v [S2\ acting on the set of vertex labels. Note that (15.361) looks slightly different 
from (I5.34p because Ei is a permutation, but E is not. The action of substituting i — > u(i) 
for some permutation u G S iv can be achieved by conjugation for Ei but not for E . 



5.5 Number of Feynman graphs from strings on a cylinder 

In this subsection we show that the formula (15.241) for counting Feynman graphs in </> 4 
theory is computing an observable in S n TFT. Recall from section [2~T1 that the expectation 
value of the observables tr(ailll),tr(a 2 U 2 ) on the cylinder with boundary holonomies U\ 
and U 2 are 

ZK<7 2 ) = / dUxdUMaxUlM^U^ZiUu^) = Y.^xl^r 1 ) (5.37) 
J 7 
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Summing this expectation value over permutations in the subgroups Hi, H 2 



Z(Hi,H 2 ) 



1 



Y Y z (vuv2) 



Hi\\H 
1 



Y Y Y 5 (^ a ^ a x ) 



(5.38) 



H\\H 



we recover the formula ( 15.211) for counting Feynman graphs in 4 theory. 

There is an interesting subtlety we should comment on. In the large N expansion of 
U(N) 2dYM, we encounter observables which can be parametrized using permutations. 
For observables constructed from gauge-invariant polynomials of degree n in the holonomy 
U, the partition functions are functions only of the conjugacy class of a in S n . In the 
above expression f)5.38j) . the sums over ^ do not run over entire S n conjugacy classes, but 
are restricted to H. 

As explained in section 2.2 the observables in the — > 1 approximation of large TV 
2dYM, for Riemann surfaces with or without boundary, can be expressed in terms of 
physical observables in S n TFT. The quantity in (15.381) is a generalization of the observ- 
ables one gets from the large N expansions of U(N) 2dYM, but it is still an observable 
in the lattice S n TFT. The boundary observables are not invariant under conjugation by 
S n elements at the boundary. The S n conjugation symmetry is broken to Hi and H 2 re- 
spectively. There is a general principle of Schur-Weyl duality which relates the symmetric 
group constructions to unitary groups (see [HI [351 ESI [37] for applications in gauge-string 
duality). The wreath products of symmetric groups have also appeared in connection 
with symmetrised traces in the context of constructing eighth-BPS operators [32] • We 
therefore expect that these more general S n TFT observables can also be expressed in 
terms of some construction with gauge theory involving unitary groups. We will leave 
this clarification for the future. 

As mentioned in Section 2, the S n TFT is closely related to covering space theory. In 
the next section we will encounter expressions of the form (I5.38P but without the sum 
over a. We will show how cutting and gluing constructions of covering spaces use the 
data jUi, H2, cr. 

5.6 The symmetry factor of a Feynman graph from strings on 
a cylinder 

We have already argued that the Feynman graphs of <j) 4 theory correspond to elements of 
the double coset, 



(5.39) 
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The S4 V on the left is the diagonal subgroup of the product group, i.e pairs of the form 
(er, a). The symmetry factor of a Feynman graph corresponding to an orbit with repre- 
sentative (<7i, 02) is the size of the stabilizer group. This can be computed by calculating 



(5.40) 



7£S 4 v m£Sv[S4 n2eS 2v [S2] 



Using one of the delta functions to perform the sum over 7 we obtain 



(5.41) 



Ati€S„[S4] [S2] 



By defining a = a 1 a 2 , we can write the formula for the symmetry factor as 



Sym(or) = Y Y 5 G"i lf W7 X ) 



(5.42) 



This expression relates more directly to the equivalent description of the double coset 



Comparing to (12.1 Op we see that this is an observable in S n TFT, with Hi,H 2 observables 
on the boundaries. This is illustrated in Fig. [T^J 

Consider a cover of the cylinder. Choose a point on one boundary circle and label the 
inverse images of that point {1,2, •• • , n}. Following the inverse images of these points 
along the boundary circle leads to a permutation fii which is constrained to be in Hi. 
Likewise there is a permutation fi 2 of {1', 2', • • • , n'} at the other boundary. Following a 
path on the cylinder which joins the two points will produce a permutation a. We are 
fixing a to lie in a fixed element of Hi \ S n /H 2 . 

From the point of view of S n lattice TFT, the line joining the two chosen points on 
the boundary circle, associated with a fixed permutation a is a Wilson line. Given the 
close relation, between S n lattice TFT and Hurwitz spaces, we can map the Wilson line 
observable to a construction in covering spaces of the cylinder. 

First notice that the cylinder can be obtained by gluing two ends of a square [AA'B'B] 
in Figure [T3J This is a topological quotient which identifies the edge AB with A'B'. To 
construct a covering space associated with the data /ii, a, fi 2 subject to fi^a^cr^ 1 = 1, we 
consider cutting the square further into the rectangles [AA'C'C] and [DD'B'B] (Figure 
[T4"l) . The cylinder is recovered by the gluings (Fig 14) 



<^M> as 



(5.43) 



AC = A'C 
DB = D'B' 
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Figure 12: The symmetry factor of a graph is an open Wilson line S n lattice TFT 
observable 
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AB = A'B' 
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Figure 13: Gluing a square to get cylinder 



CC = DD' (5.44) 

Now take n copies of these pairs of rectangles, with labels [AiA^C-Ci] and [DiD^B^B^, 
with % ranging from 1 to n, as in Figure [15j For the application to (I5.38P we have n = 4v. 

The union of rectangles {[Ai^C^Ci], • ■ ■ , [A n A' n C' n C n }} is quotiented by identifying 
edges AiCi to the edges A\C[ by the permutation fix. 

(5.45) 

The rectangles [DiD^B^Bi] are quotiented by identifications 

D 'i B 'i = D^B^ (5.46) 

Finally we identify 

D % D i = awCw (5-47) 
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Figure 14: Gluing a pair of rectangles to get cylinder 



The covering map is specified by mapping each of the n rectangles, to the original rect- 
angles without labels in the obvious way (see Figure [T5|) 

The condition /i^V/^ct -1 = 1 ensures that if we consider the inverse image of the 
closed contractible path shown in blue in Figure [161 the inverse image is also contractible. 

This type of construction is a standard part of covering space theory in connection 
with the Riemann existence theorem (see for example |39j). 

The generalizations of ( 12. 8ft visible in (15. 38ft and (15.421) involve restricting the sums 
over fix, fi 2 to specific subgroups Hi, H 2 and leaving the a unsummed. From the point of 
view of S n TFT these constraints are easily implemented, since the degrees of freedom on 
each edge form the whole group, where elements and subgroups can be chosen. However 
the boundary observables are not invariant under conjugation by S n and have no dual in 
the large N 2dYM that we yet know how to construct. Given the close connection between 
S n elements and monodromies of covers, it is not surprising that a definite covering space 
construction for specific /ii,/i2,cr exists. We have provided such an explicit cutting-and- 
gluing construction here. There are clear analogies between the construction given and 
the implementation of twisted boundary conditions in D-branes on orbifolds [3D]. The 
presence of a Wilson line defect located at a line joining the 2 boundaries of the cylinder 
in the S n TFT picture suggests that there should be a D-brane interpretation. We leave 
this as a future direction of research. Developing the connection of the lattice S n TFT to 
the axiomatic approach to TFT and branes in [36] maybe a useful approach. 

In the standard discussions of the Riemann existence theorem for (branched) covers, 
one has a two-way relation. On the geometrical side, there are maps / with the equiva- 
lence / = / o (p for holomorphic automorphisms <fi (see more details in P, [25]). On the 
combinatoric side there are permutations with conjugation equivalence. From permuta- 
tions we can construct covers, and vice versa. The equivalence classes map to each other. 
Here we have permutations, but with some restrictions having to do with choices Hi,H 2 . 
We can still construct some covers with this data. We have not fully articulated the cor- 
rect equivalences on the geometrical side for a 1-1 correspondence. The construction we 
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Figure 15: Gluing copies of rectangle pairs to get covering of cylinder - identifications 
determined by fii, fi 2 , & 

gave should provide useful hints for the precise definitions on the geometrical side, which 
will provide the equivalence. We leave this as a problem for the future. 

6 The number of Feynman graphs 

In section 15.31 a formula expressing the number of Feynman graphs in terms of the star 
products of two cycle indices was given. This result is useful because formulas for the 
cycle index of the wreath product of two symmetric groups are known. Using these results 
we will write down rather explicit formulas for the number of vacuum graphs in 4 theory 
in section 16.11 
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C — c 
D D' 

B B' 

Figure 16: Closed contractible path 



6.1 Some analytic expressions exploiting generating functions 

Let di denote the total number of vacuum Feynman graphs in 4 theory with % vertices. 
To obtain an analytic formula for di we will need the cycle indices of two wreath products, 
SnfSy and S n [S2\- The known generating functions for these wreath products are 



Z s ^[t,x] = J2t n Z Sn[S2] (x) = e ££i (6.2) 

n 

To compute N(Z(Hi) * Z(H 2 )) take the product of the coefficients of Yii % Pi in the two 
cycle indices, and then further take the product with rL^Pi- This can be accomplished 
if we make the substitution %i — > y/ii/i, to obtain 

Z s °°^ [t, yi ] = Z s ^ [t, Xi = Vt Vl ] 

Z s °° [Si] [t, Vi] = Z s °°^ [t, Xi = ViVi] (6.3) 

and then replace the {yii)i) Pi with pi\. This leads to the following formula for the number 
of 4> A vacuum graphs with v vertices 

d v = f^flj ^e-^^^^^,^]^^^,^] (6.4) 

The integrals over i/i and tji ensure that only terms with equal powers of yi and i/i con- 
tribute, and they implement the substitution (yiy~i) Pi — > Pi\. The contours of integration 
over ti and t 2 are both counter clockwise and they both encircle the origin. 

We could contemplate many ways of refining this result. For example, can we deter- 
mine how many of these vacuum graphs are connected? Before turning to this question, 



32 



it is instructive to ask how we can identify if a graph is connected or disconnected from 
its description in terms of the pair S , S x . A Feynman graph (S , £1) with v vertices is 
disconnected if for some subgroup G = S^-p x S p with p > we have cr(Si) G G for all 
a G Aut(So). Let q denote the number of connected Feynman graphs with z vertices. 
The total number of vacuum graphs with % vertices is the coefficient of g % in the partition 
function 

oo oo 

2=1 1=1 

This formula, given for example in [20] , can be used to determine the q once the di have 
been computed. It is now straight forward to determine the number of vacuum diagrams. 
See Appendix [D] for numerical results. 



1-g* 



(6.5) 



6.2 Scalar fields with external legs 

The next natural generalization is to consider Feynman graphs with E external legs. 
Summing Feynman graphs with E external legs produces an £"-point correlation function. 
For the generic case when all E'-points corresponds to different spacetime events, graphs 
obtained by permuting labels of the external points give distinct contributions to the 
correlation function. For this reason the automorphism group of the graph does not 
include any elements that permute the labels of external legs. 

For a graph with v vertices and E external legs, we obtain a total of 4v + E half 
edges. Ei is now an element of the conjugacy class [2 2 " + 2 £: ] of S iv+ E, while S contains v 
4-tuples of numbers specifying how the half edges connecting to the vertices are labeled, 
as well as E numbers specifying how the half edges connecting to external points are 
labeled. Since we do not have automorphisms that permute the half edges connecting 
to external points, Hi = Aut(S ) = S^Sy x (Si) E . We can also write the subgroup as 
iSyfiSi], with the understanding that it is acts by keeping fixed the last E integers from the 
set {1,2, •• • ,4v + E} that S iv+ E acts on. Consequently, Feynman graphs in 4 theory 
with v vertices and E external points are in one-to-one correspondence with elements of 
the double coset 

{S V [S 4 ] x Sf ) \ S 4v+E /S^ +E) \S 2 ] (6.6) 

2 

To compute the number of graphs the only new cycle index we need is 

z *[*]xsf (xi) ... )Xn) =x f Z *[*]( Xl> ... jXn ) (6.7) 

It is now straight forward to obtain the number of Feynman graphs in 4 theory with E 
external lines 

Af(Hi, H 2 ; G) = N{Z{Hi) * Z(H 2 ); S iv+E ) = £ Z^* s ? Z^^Symip) (6.8) 

p\-(iv+E) 
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The notation p h (4v + E) indicates that p runs over all partitions of 4v + E. 

We can again refine this counting by asking how many of these graphs are connected. 
For concreteness, focus on the case E = 2. Let d 2t i denote the number of Feynman graphs 
with two external lines and % vertices and let c 2j i denote the number of these graphs that 
are connected. We can again write down a partition function which gives the total number 
of graphs with two external lines and % vertices as the coefficient of g l . In this formula the 
number of vacuum graphs q and di defined in the last subsection participate. We find 



n 

i=X 



C2 > kgk ) = [Yl d2 > igi 



(6.9) 



,fe=0 



,fc=i 



This formula can be used to determine the c 2 j once the d 2 ^ have been computed using 
( 16. 8p . See Appendix |D1 for numerical results. 



6.3 Generalization to (ff theory and other interactions 

Although our discussion has focused on 4 theory it should be clear that our methods are 
general. To illustrate this we will now consider a <fr 3 interaction. A graph with v vertices 

3v 

has 3v half edges. Si is now an element of the conjugacy class [2~] of Ss v , while S 
contains v 3-tuples of numbers specifying how the half edges connecting to the vertices 
are labeled. We have Hi = S V [S3] and H 2 = Sav [S-j] ■ Vacuum Feynman graphs in 3 
theory with v vertices are in one-to-one correspondence with elements of the double coset 

S v [S 3 ]\Sto/Si»[S 2 ] (6.10) 

To compute the number of graphs the only new cycle index we need can be read from the 
generating function 

oo 

Z s °°^ [t; Xl , x 2 , ■ ■ ■ ] =J2 t n Z s -^ ( Xl ,x 2 , ■ ■ ■ , x n ) 

n=0 

= e E£iw( 3, ?+ 3!B M**+ 2 *si) (6.11) 
It is now straight forward to obtain the number of vacuum Feynman graphs in 3 

^3 [$2] 

N(Z(H 1 )*Z(H 2 )) = J2 Z p v[S3]z p ¥ 'SymW (6-12) 

p\-3v 

Feynman graphs with E external legs are in one-to-one correspondence with elements 
of the double coset 

(S V [S 3 ] x Sf ) \ S 3v+E /S {3v+E) [S 2 ] (6.13) 
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Using the generating function 

Z Sn ^ xS ?{x u --- ,x n ) = xfZ Sn[S3] (x u --- ,x n ) (6.14) 

it is straight forward to compute the number of Feynman graphs with v vertices and E 
edges 

^3 E [^2] 

Af(H 1 ,H 2 ;S 3v+E )) = N(Z(H 1 )*Z(H 2 );S 3v+E )= ^ Z $»W xS ? Zp^ °~ Sym(p) 

ph(3v+E) 

(6.15) 

The notation p h (3v + E) indicates that p runs over partitions of (3v + E). See Appendix 
ID1 for numerical results. Cubic graphs have recently played a role in studies of N = 8 
SUGRA [J9j as well as in the classification of N = 2 4-dimensional gauge theories |48j . 
The 3 theory in 6 dimensions is known to be asymptotically free [5U]. Hence this case is 
of special interest. 

Finally, consider a theory with both cubic and quartic vertices. To count the vacuum 
graphs having v 3 cubic and u 4 quartic vertices, we would repeat the above discussion with 
S V3 [S 3 ] x S V4 [S 4 ] replacing the groups S V3 [S 3 ] (for <p 3 ) or S V4 [S A ] (for </> 4 ). 

6.4 Interpretation in terms of 2dYM string 

We have already argued that the number of vacuum graphs in 4 theory is computing an 
observable in the Lattice TFT with S n gauge group. The counting of Feynman graphs 
in this section also has a formulation in Lattice TFT: One simply uses (I5.38P with the 
appropriate Hi,H 2 described above. 

The formula (15.421) for the symmetry factor as an amplitude is also directly applicable 
here. 

7 Feynman graphs as orbits: a view of ribbon graphs 

The identification of Feynman graphs as elements of the double coset S n \ (S n x S n ) / (Hi x 
H 2 ) is a nice unifying picture. The groups Hi and H 2 are the symmetries of the interaction 
(black) and bivalent (white) vertices respectively. This allows a simple generalization from 
ordinary graphs to ribbon graphs: symmetries must by definition preserve the genus of 
the ribbon graph. The cyclic order of the labels at a black vertex must be respected if the 
genus is to be preserved. Thus, for example, the replacement Hi = S^SV] — > Hi — S^Zzi] 
takes us from counting Feynman graphs in <fr A theory to counting ribbon graphs in Tr0 4 
theory. 
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Feynman Graph Vertex 
S 4 symmetry 



Ribbon Graph Vertex 
Z 4 symmetry 



Figure 17: Cyclic versus symmetric for Feynman vertex versus ribbon vertex 



It is now clear that ribbon graphs can also identified with elements of a double coset. 
For example, the vacuum graphs with v vertices of a matrix model with Tr0 4 interaction 
are elements of the coset 



where we have identified H\ = S v [7,i] and Hi = S^fSV]- 

7.1 Number of ribbon graphs using commuting pairs 

Using the insights of the previous subsection and repeating the argument of section El 
we can generate all possible ribbon graphs by allowing S I „[Z 4 ] to act on the set of all 
possible Wick contractions, that is, the conjugacy class [2 2v ] of S iv . Each orbit of S , „[Z 4 ] 
is a distinct ribbon graph. Thus, the number of ribbon graphs can be obtained, using 
Burnside's Lemma, as 



(7.1) 



Number of ribbon graphs 




1 



E 



(7.2) 



Hi 



a£S iv /H 2 



Using the results of section 15.41 we can also write 



Number of ribbon graphs 




7G52„[S 2 ]ctG[4"] 




(7.3) 
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There is also a way of writing this that is manifestly symmetric between exchange of 
Hi, H 2 

Number of ribbon graphs = 1 V V" V" Xj?(7i)X/?(72) 

LHi iz2 n! z — ' ^ — ' 

= £^£^5* ( 7 - 4 ) 

Recall the notation Adf for the multiplicity of the one-dimensional representation of H\ 
when the irrep R of S± v is decomposed into representations of the subgroup H±. 



7.2 How many ribbon graphs for the same Feyman graph ? 

Ribbon graphs come equipped with a natural notion of a genus. For this reason not 
all of the contractions leading to a given Feynman graph will produce the same ribbon 
graph. As an example, consider vacuum graphs in c/> 4 theory with v — 1 vertex. All three 
contractions give the same Feynman graph, but two of them give a genus zero ribbon 
graph and one gives a genus one ribbon graph : see Figure [18j 

The line of thinking developed above allows an elegant answer to the question of how 
many ribbon graphs are there for a given Feynman graph. In this case So can be chosen 
as a permutation made of 4-cycles 

S = (1234)(5678)...(4u-3,4i;-2,4i;- l,4v) (7.5) 

The symmetry group of this permutation S^Z^. This is the set of permutations in S& v 
which commutes with So- The Wick contractions are described by Si which are pairings, 
i.e permutations in the conjugacy class [2 2v ] of S^. Ribbon graphs are orbits of S^fZj 
action (by conjugation) on the permutations in [2 2v ]. We can first decompose the set 
of all possible the pairings in [2 2v ] into orbits of the group S^SV]. Each such orbit is 
a Feynman graph. Then decompose each orbit into orbits of the subgroup S^Z^ of 
S v [S4] . The elements within one orbit of S v [S4] will fall generically into multiple orbits of 
S^fZj, corresponding to multiple ribbon graphs. We can view the genus of the worldsheet 
described by the ribbon graph as an invariant of the orbits of S , „[Z 4 ] acting on 2-cycles. If 
two 2-cycles are in the same orbit, they must be associated with the same genus. Consider 
the permutation obtained by multiplying So with any Wick contraction which belongs to 
a fixed orbit of S^fZ/t]. Denote this permutation by cr 3 and the number of cycles in this 
permutation by c. The genus g of each S^fZ/t] orbit is related to the number of cycles in 
the permutation c and number of vertices v by 

2-2g = c-v (7.6) 

This is obtained by applying the Riemann Hurwitz formula to maps the case of maps 
with three branch points determined by S , S 1; 03 |19j . 
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It is straightforward to perform, in software such as GAP, the decomposition into 
orbits of iSjiSy, and then refine the decomposition according to S , „[Z 4 ]. 



3 



> 2 




+ 




Figure 18: At v — 1 the three possible Wick contractions give the same Feynman graph. 
Two of them give a genus zero ribbon graph and one gives a genus one ribbon graph. 



A number of interesting directions can be contemplated. Different ribbon graphs 
with the same underlying Feynman graph will have the same space-time integrals but 
possibly different genus. Our methods provide information about how the space-time 
dependence allows a certain range of genera. When we are dealing with vacuum graphs, 
each Feynman graph contributes a number. However, our techniques can be applied to 
graphs that have external legs in which case we have explicit space-time dependences. 
Presumably a Feynman graph with a certain "complexity" - reflected in its space-time 
dependence, will allow a certain range of genera. The more complex it becomes, the more 
genera it will allow. This could be studied quantitatively with the current set-up. 



7.3 Using cycle indices for ribbon graphs 

In this section we would like to count the number of ribbon graphs using cycle indices. 
For 4-valent ribbon graphs, the total number of graphs for v vertices, is obtained by using 
the formula 

N( Z(S v [Z 4 ])*Z(S 2v [S 2 }) ) (7.7) 
Thus, we need the cycle index of S , l ,[Z 4 ]. The cycle index of any cyclic group is 



Z £n (x) 



n 



where <p(d) is the Euler totient function. For Z 4 we find 

1 



Z M (x) = j(xj + x 2 2 + 2x 4 ) (7.9) 
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Using known results for the cycle index of wreath products we now find 




(7.10) 



n 



It is now straight forward to obtain explicit answers for the number of ribbon graphs. 

The methods at hand also apply to ribbon graphs with arbitrary numbers and types 
of traces. For a ribbon graph with v\ single traces, v<i double traces, v% triple traces etc. 
we would use the group S Vl x S V2 [Z 2 ] x S V3 [Z 3 ] • • • to count graphs. 

8 QED or Yukawa theory 

We have focused on real fields corresponding to particles which are neutral. In this section 
we explain how our methods can be extended to the case of complex fields corresponding 
to charged particles. We will refer to the charged particle as an electron with a view to 
applications to QED. Since are just counting Feynman graphs, we do not track minus 
signs coming from the anti-commuting nature of fermion fields. As far as the number of 
diagrams or their symmetry factors goes, there is no difference between QED or Yukawa 
theory. Of course, in QED certain diagrams vanish automatically due to Furry's Theorem. 
In the next section we will consider the problem of counting the QED diagrams that remain 
after Furry's Theorem is applied. 

The fact that we are dealing with charged particles implies that each edge will now 
have a preferred direction, determined by the flow of charge. The model we have in mind 
has a cubic vertex in which the charged particle interacts with a neutral particle. This 
structure of the vertex matches both the Yukawa interaction and the coupling of charged 
fermions to a gauge field. We will refer to the neutral particle as the photon. Thus, each 
vertex has 3 edges: an incoming electron, an outgoing electron and a photon. Consider 
vacuum Feynman graphs made of 2v vertices. We again clean the graph, by introducing 
new vertices in the middle of each edge. Label the vertices so that 1, 2v are attached 
to the incoming electrons, 2v + 1, • ■ • , 4w go with the outgoing electrons and 4v + 1 ■ ■ - 6v 
go with the photons. 

So we can describe the vertices by 



2r 



(8.1) 



i=l 



The Wick contractions are pairings of the form 



Si = ]J{i,2v + -a x 



(8.2) 



i=l 
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Figure 19: QED vertices 



The permutation <j\ is in [2 V ] inside the S 2v which acts on the last 2v among the {1, • • • , 6v}. 
The permutation T\ is a general permutation of 2v objects. It tells us which of the first 
2v (incoming electrons) go with which of the 2v + 1 • • • Av (outgoing electrons) . The sym- 
metry 7 of S is just S 2v of exchanging the 2v brackets. Because all legs are distinct there 
is no symmetry exchanging the legs at a vertex. Given the form of S we see that this is 
the diagonal S 2v of the S 2v x S 2v x S 2v acting on the electron, anti-electron and photon 
labels. We will call this Diag 3 (S , 2„). 

The automorphisms of the Feynman graph are those elements 7 G Aut(S ) = Diagg^^ 
which also leave fixed the Si pairing determined by (a\, 7~i). The order of this group again 
determines the symmetry factor of the Feynman graph. 

Distinct Feynman graphs are orbits of the permutations 



Diag 3 (S 2v ) -> S 2v x S 2v x S 2v 



;.3) 



acting on Si(a"i, 77), where the a±, T\ are described above. We know that the stabilizer of 
o\ G [2 V ] in S 2v is conjugate to S V [S<2]. This leads to 



S 2v /S v [S 2 ] = [2 V ] 



5.4) 



The permutation n mixes the first 2v incoming electrons with the next 2v outgoing 
electrons. The stabilizer of such a permutation is the diagonal S 2v which simultaneously 
moves the first 2v with the next 2v, so that T\ is running over the coset 



(S 2v x S 2v ) /Di&g 2 (S 2v ] 



!.5) 



We use the subscript 2 because it is the diagonal of the two S 2v . This is enough to 
conclude that the Feynman graphs are in one-to-one correspondence with the points of 
the double coset 



Di&g 3 (S 2v ) \ (S 2v x S 2v x S 2v ) I (Di&g 2 (S 2v ) x S V [S 2 ]) 
The Diag 2 is the diagonal of the first two S 2v in S 2v x S 2v x S 2v . 



1.6) 
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8.1 Double cosets of product symmetric groups 

To count the number of Feynman graphs, we need to count the number of points in the 
double coset 

Hx \ (S 2v x S 2v x S 2v )/H 2 (8.7) 

where Hi = Diag 3 (S 2v ) and H 2 = Dia.g 2 (S 2v ) x S V [S 2 ]. It is instructive to approach this 
problem by first considering cosets 

#1 \ G/H 2 (8.8) 

where G is a general product of symmetric groups S ni x S n2 • • • x S nk , and Hi, H 2 are 
general subgroups. Let us denote by 

M{H U H 2 -G) (8.9) 

the number of points in this double coset. The case J\f(Hi, H 2 ; G = S n ) = N(Hi * H 2 ) is 
the one we already discussed in sections 15.21 and 15.31 In the case at hand, we need 

M(H u H 2 ;S 2v xS 2v xS 2v ) (8.10) 

As a consequence of the fact that the numerator group has changed, when we compute 
the star product, it will turn out that we should multiply the coefficients of Z(Hx), Z(H 2 ), 
and now weight them by the symmetry factors appropriate for S ni x S n2 x ■ • • S nk . This 
is a generalization of Read's original formula [20] so it is worth discussing in some detail. 

When the numerator group G is a product of symmetric groups, a subgroup such as 
Hi will have some number of elements in each conjugacy class of G, which is specified by 
three partitions. 

{ p W jP W ... p ( fe )} (8.11) 
For the case at hand k = 3. The symmetry of such a conjugacy class is 

3 2v 

Sym({p«}) = J] Sym(p«) = n \[f?pf\ (8.12) 

i i=l j=l 

We will restrict the discussion to k — 3, but the generalization to any k is immediate. 
For a subgroup H of G it is natural to define 

Z H -°= E ^^ll^yf'zf (8.13) 

p(l),p(2) )P (3) i,j,k 
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where the coefficients Z ^ 2 ) p <3) keep track of the number of permutations in the sub- 
group with the cycle structure specified by the 3 partitions. For two subgroups Hi,H 2 , 
we can use the two cycle indices (with respect to G) to define a generalization of the star 
product in ( 15.1 9p as follows 

M{H X , H 2 ; G) = J2 ^S >p (3)^S jP (3)Sym({p«}) (8.14) 
p 

We can understand this formula by adapting the reasoning in [20] . For any p G G, the 
product hip h 2 gives \Hi\ x \H 2 \ elements, that all belong to the same equivalence class in 
the coset (I8.9p . as hi runs over Hi and h 2 runs over H 2 . Some elements will be repeated. 
When this happens 

hiph 2 = hiph 2 (8.15) 



which implies that 

h^hi = ph 2 h 2 x p^ (8.16) 

Clearly hi hi G Hi n pH 2 p~~ x . Given hi, h 2 and an element of Hi D pH 2 p~ 1 , hi and h 2 
are uniquely determined. Consequently the number of times that p is repeated, v(p), is 
\Hi n pH 2 p~ 1 \. Every element in G which appears in the equivalence class of p comes 
with this same multiplicity in HipH 2 . Hence, the number of distinct elements in the 
equivalence class of p 

, . H\ Hi , _. 

*» = -iw 1 (8 ' 17) 

Consider 

J2"(P) (8-18) 

For an equivalence class containing p, we get a contribution n(p)u(p) = \Hi\\H 2 \. This is 
independent of the equivalence class. If there are N d equivalence classes, we get 

J2"(P) = N d \Hi\\H 2 \ (8.19) 

p<=G 



Rearranging we find 



N d = ■ 

Hi H 2 



LE^p) (8-20) 



peG 
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To compute the sum over p of v{p) we can choose U\ G H\ and u 2 G H 2 and count how 
many p's obey 

mi = pu 2 p~ X (8.21) 

Since iti and u 2 are conjugate, they have the same cycle structure. Thus, to compute the 
sum over p of v(p) we need to fix a cycle structure, multiply the number of elements in 
Hi with this cycle structure by the number of elements in H 2 with this cycle structure, 
then multiply by the number of elements of G that leave this cycle structure invariant, 
and then finally sum over cycle structure. This is precisely what (I8.14p is computing. 



8.2 Analytic expressions for number of Feynman graphs in 
QED 

In this section we will apply the formulas obtained in the previous section to obtain 
explicit results for the number of Feynman graphs. We need two cycle indices. For 
Hi = Diag 3 (S' 2 „) we have 

1 2v 

ph2v J v ^ y j=l 

For H 2 = Diag 2 (S , 2l) ) x S V [S 2 ], we have 

z „^a = J2J2 * 2*W n^)^ (8.23) 

qh2vr\-2v y W j=l 

To calculate the J\f(Hi, H 2 ;G), we need to multiply like terms in the f )8.22p and ( I8.23p . 
weighted by an appropriate symmetry factor. Picking up like terms forces q = p = r so 
that 

i r^S v [S2] 



Z^ [S2] Sjm{p) (8.24) 

p\-2v 



For the wreath product 

CO 

Z s °°^[t;xi,x 2 ,---} = J2 ZSn[S2] ( x i>--- > x *n)t n 



n=0 

e J2T=i ( 8 .25) 
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To get the desired counting from this, we need to replace Yli X T w ^h i Pi Pi\. Equivalently, 
do a replacement X{ —> iy^ and expand in replacing y\ l with pj.. In terms of 

Z[t- Vl ] = Z s °°lt;xi = i yi ] 

= 6 J2ZitM+^) (8.26) 

the counting of vacuum graphs with 2v vertices can be written as 

—r(\ll d yi e^)Z[t ]yi ,y 2 ---} (8.27) 
See Appendix [D] for explicit numerical results obtained using these formulas. 



8.3 QED and TFT for product symmetric groups 

The number of Feynman graphs is counted by the formula 

jjfj^ E E E^wr" 1 ) (8-28) 

We have already established that 

G = x S 2v x 5*2^ 
if x = Diag 3 (S 2v ) 

H 2 = Diag 2 (5 2 „) x S V [S 2 ] (8.29) 



This is a cylinder partition function in x x gauge theory which is Schur- 
Weyl dual to the large N limit of U(N) x U(N) x U(N) gauge theory. As before the 
correspondence holds in the zero area limit. 



9 QED with the Furry's theorem constraint 

Above we have counted the total number of Feynman diagrams in QED or Yukawa theory. 
In the case of QED, Furry's Theorem proves a fermion loop punctuated with an odd 
number of photons vanishes [H]. In this the number of QED Feynman graphs, not 
vanishing by the Furry constraint, are counted. Enumeration of some Feynman graphs 
with their symmetry factors for this case is given in [43J. 
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T=(12) : non-zero 
by Furry's Theorem 



T=(1)(2) : zero by 
Furry's Theorem 



Figure 20: Even cycle lengths and Furry's theorem 

9.1 Furry's theorem and a constraint to even cycles 

Following the discussion of the last section, each QED Feynman graph is specified by 

So = [ < i; 2v + i; 4v + % > 

i=i 

Ziiaun) = JJ(z,2 v + n(i)) • U\ (9.1) 
i=i 

Recall that T\ specifies how incoming and outgoing electrons are connected. Consequently, 
the number of cycles appearing in n is equal to the number of fermion loops in the graph 
and the length of each cycles is equal to the number of photons decorating the loop. It 
is clear that the constraint that each fermion loop has only an even number of photons 
is easily imposed in the (E , Ei) description by requiring that t\ consists of permutations 
which only have cycles of even length, as shown in Fig 20. 

The orbits of the group G = S 2v x S 2v x S 2v acting on the pairs (E , Ex) correspond 
to the inequivalent Feynman graphs. The set of 7 G G which preserve the pair (Eo, Si) 
generates the automorphism group of the graph and the order of this automorphism group 
is the symmetry factor. Preserving Eo forces 7 to be in the diagonal of G : 

For 1 < % < 2v 
i — > 7(2) 
2v + % -> 2v + 7(z) 

4v + i 4w + 7(z) (9.2) 
We can write this as 70707, where the first 7 acts on the first 2v according to the second 
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line of (19. 2p . while leaving the subset {2v + 1, • • ■ , 6v} fixed ; the second 7 moves only the 
subset {2v + 1, • • • , 4t>} according to the 3rd line; and the third 7 moves only the subset 
{4v + !,-■■ , Qv} according to the 4th line. 

The symmetry factor of a Feynman graph specified by the standard E and a general 
Ei is 

|Aut([Ei] Hi )| = ^ 5(7 07o 7 E 1 7- 1 o7- 1 o7- 1 Sr 1 ) (9.3) 

yeS 2v 

Ei is not a permutation in G, since the T\ mixes the first 2v with the second 2v. The 
multiplication can be viewed in Sq v or Sa v x S 2v 

A simple application of the Burnside Lemma implies that the number of Feynman 
graphs remaining after Furry's Theorem is applied is 

T^ry Y Y Y 5( 7°7°7 EiOi, r 1 ) 7 ~ 1 o o 7 ~ x E^Vi, n) 

^ >' -yeSzv <n€[2 v ]es 2 v riG5 2t ,:[Ti] even 



0- 



Another approach to the same counting problem would be to implement the Furry 
constraint in the double coset language. Recall the double coset relevant for QED is 

Diag 3 (5 2 „) \ (S 2v x S 2v x S 2v )/Dmg 2 (S 2v ) x S V [S 2 ] (9.5) 

Using (I5.2ip with the G,Hi,H 2 identified according to (19. 5 h the size F of this double 
coset is counted by 

(2v)\(2v)\2 v v\ 

E E E E 

v rK r ' aeS 2v t£S 2v p€S v [S 2 ]h,b 2 ,b 3 eS 2v 
$S 2v xS 2v xS 2v (c o- o- hob 2 ob 3 tot op b^ob^ob^ 1 ) (9.6) 

b\ acts on the incoming electrons, b 2 on the outgoing electrons and 63 on the photons, a 
permutes vertices in the graph, r permutes electron lines and p permutes photon lines. 
This can be simplified to 



E E E E 



(2v)\(2v)\2 v v\ 

v ' v ' veS 2 v t€S 2v peS v [S2]bi,b2,b 3 eS 2v 
Ss 2v xs 2v ( a a h°b 2 tot b^ob? 1 ) 

Ss^hpbz 1 ) (9.7) 

Solving the second delta function for a, we can plug back into the first delta function and 
do a redefinition of the sums Y^ fei b2 to absorb the 63, so as to get 

f = - y v v 

(2v)\2 v v\ ^ ^ ^ 

y ' tgs 2v P eS v [S 2 ]bi,b 2 eS 2v 
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$s 2v xS 2v (p°P b 1 ob 2 tot b^ob^) 
= (2vY2"vl £ £ £ S s 2v (pbirb^)5s 2v (pb2rb^) 

= ^ £ £ wWp- 1 ) (9.8) 

In the last line we have recognized that, given the interpretation of 61,62 it is clear that 
Ti = b~ 1 b 2 is the permutation t\ in (19. ip . The Furry constraint is now easily implemented. 
Thus, the number of Feynman graphs remaining after the Furry constraint is implemented 
is given by 

Number of Feynman graphs for QED with Furry constraint 
= 2^i £ £ WW," 1 ) (9.9) 

ri£S 2v :even pGS v [S 2 ] 

The permutation 73. is constrained to have even cycles only and we are summing over 
elements of S V [S2\- For each element of 5^ [£2], the weight is the number of permutations 
with even cycles only, which commute with the given permutation in S^fS^]. 

From the first line of ( 19. 8p we see that an equivalent double coset description of the 
counting is 

S 2v \ (S 2v x ^VDiag^^]) (9.10) 

This gives a slightly simpler connection to observables for 2d Yang-Mills with cylinder 
target space. Namely we have a connection to U (00) x U (00) rather than U (00) x U (00) x 
U (00) as described earlier. 



9.2 From QED to ribbon graphs 

Comparing the expressions (19.81) and (19.91) with (17.41) for ribbon graph counting, it becomes 
clear that the counting of QED Feynman graphs can be matched with counting ribbon 
graphs. The restriction to vertices of valency 4 in (17.41) is being relaxed to allow arbitrary 
valencies in (I9.8P and arbitrary even valencies in (19.91) . We will now explain the bijection 
between ribbon graphs and QED/ Yukawa graphs which explains the equality of counting 
formulae. 

The general QED /Yukawa vacuum graph has loops with vertices where the photon 
joins the loop. These loops determine 2v cycles which form a permutation in S^. To 
build a bijection to ribbon graphs think of the photon labels in Figure [19] as attached 
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to the vertices. The photon contractions determine an element of [2 2v ]. At the centre 
of each electron loop draw a point and radiate edges (spokes) to intersect the edges of 
the loop, one spoke between each pair of electron-photon vertices. Use the arrow on the 
electron propagator to move each end of the photon propagator along the electron loop 
towards the intersection of the spoke. Erase the electron loop, leaving the spokes and 
the photon propagators connecting them. This gives a graph with vertices of arbitrary 
valency, but with the vertices being equipped with a cyclic order, which is nothing but a 
ribbon graph, completing the construction of the bijection. This process is illustrated in 
Figure [2TJ The bijection demonstrates that the number of QED vacuum graphs with 2v 
vertices is the same as the number of ribbon graphs with v edges (the photon propagators 
become edges). 




Figure 21: QED graphs to ribbon graphs 



This means that the total valency of the vertices in the ribbon graph is 2v. Each 
vertex can have any integer valency compatible with this constraint. Thus, following the 
discussion (17.31) and (17.71) . the number of QED/Yukawa vacuum graphs with 2v vertices 
is 

Number of Feynman graphs for QED = N{Z{H l ) * Z(H P )) (9.11) 

p\-2v 

where p is a partition of 2v, the group Hi is S 1 ,,^] and the group H p is 

SpJZi] x S P2 [Z 2 ] x ■ ■ ■ S P2v [Z 2v ] (9.12) 

The function N of two cycle index polynomials is as defined in (15.191) . For QED with the 
Furry constraint implemented, the valencies are constrained to be even so that 

Number of Feynman graphs for QED with Furry constraint 
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N(Z(H 1 )*Z(H p )) 



(9.13) 



ph2v 
p even 



Given the connection between ribbon graphs and Belyi pairs [19] we see that QED 
counts the number of clean Belyi pairs (bi-partite embedded graphs with bivalent white 
vertices) with degree v and black vertices of any valency. QED with the Furry constraint 
counts the number of such Belyi pairs with black vertices of even order only. 

10 Discussion 

10.1 Permutations and Strings : a new perspective on gauge- 
string dualities. 

We have seen that the counting of Feynman diagrams can be interpreted in terms of string 
amplitudes with a cylinder target space. The nature of the interaction vertices determines 
a permutation symmetry group H\ and the Wick contractions determine a permutation 
symmetry group if 2- The counting could also be written in terms of commuting pairs of 
permutations, where one permutation lives in one of these subgroups, and another one 
lives in a coset involving the other group. This latter realization can be interpreted in 
terms of covering maps of a torus target space, such as those that arise in the string 
theory of 2d Yang Mills for torus target space [3j H], except that there is a constraint on 
the two permutations associated with the cover, by following the inverse image of paths 
along the a and b cycle. 

In 2d Yang Mills, one encounters a sum over commuting permutations si, s 2 which is 
equivalent to counting the set of all covers of the torus by a torus. This sum is invariant 
under the SL(2, Z) of the space-time torus, which acts on these permutations by 



In the case at hand, the SL(2, Z) of space-time is broken in that si is summed over Hi 
while s 2 is summed over S iv /H 2 = [2 2v ]. 





and obeys the SL(2, Z) relations : 




(10.2) 
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While making two choices of observable at each end of a cylinder seems like a natural 
construction for string theory with a cylinder, the constraints on the covering maps of 
the torus look rather intricate, and are not of any sort that has been discussed in the 
literature on topological strings. Nevertheless, there might well be a string construction 
which sums over maps with constraints of this sort. 

In any case, the Feynman diagram combinatorics suggests two emergent dimensions 
of cylinder or torus, with one formulation possibly more appealing than the other. In the 
construction of correlators of general trace operators in the Hermitian Matrix model, there 
is another 2-dimensional target space emerging, which is the sphere with three punctures 

In the application of Feynman graphs to QFT, the Feynman graph counting is only 
part of the answer. When there are external edges, each edge is associated with a space- 
time position (or momentum), each of which takes values in M 4 . The computation of 
Greens functions (4>(xi)4>(x 2 ) • • -4>(xe)) is a sum over Feynman diagrams, which live in 

DC{v, E) = {S V [S 4 } x Sf ) \ S 4v+E /S 2v+E/2 [S 2 ] (10.3) 

This has an action of Se- Projecting down the Se orbits gives another double coset. 
There is a fibration 

(^[Sy x S l ) \ S^v+e/ S2V+E/2 [S 2 ] <- S E 
I 

(S V [S 4 } x S E ) \ S 4v+E /S 2v+E/2 [S 2 } (10.4) 

We may write 

G{x\, x 2 , ■ ■ ■ , x E ) = 7^ G(x!,x 2 , ■■■x E ; S) (10.5) 

E 

where £ lives in the double coset (110.31) . The Bose symmetry of G(xi,x 2 , ■ ■ ■ ,xe) arises 
as an Se invariance of the summands 

G(x ail) ,x a{2) , ■ • -x^e)] o-(E)) (10.6) 

This shows that in some sense the space-time coordinates (IR 4 )- B and the double coset 
DC(v,E) associated with covers of a cylinder may usefully be viewed as being on an 
equal footing when we work out the physics. It is tempting to infer that this is indicative 
of an underlying M 4 x Cylinder or IR 4 x Torus in four dimensional QFTs. We have shown 
that the extra two dimensions emerge from the Feynman graph combinatorics, but it 
remains to be seen if there is a dynamical six dimensional picture along these lines, and 
how it might relate to other appearances of six dimensions in 4D QFTs such as [5T~] . 

A very interesting problem is to interpret the full correlators or S-matrices, in terms 
of the cylinder (or torus) emerging from combinatorics, with the R 4 which is manifestly 
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there to begin with. This would be an example of a holographic dual of a quantum field 
theory, analogous to AdS/CFT, where from the point of view of the dual gauge theory, 
an emergent S 1 along with an emergent radial direction, resulting in M 5 x S 1 , with the 
original R 4 of the gauge theory realized as the boundary of the M 5 . 

If this is indeed possible with 4 theory, it would be an example of gauge-string duality 
with an interesting difference from all the examples known so far. It would not rely on 
the large N expansion, and ribbon graphs thickening into string worldsheets. Rather it 
would be a realization through the fundamental link between permutations and strings, 
of the physical expectation of holography [521 [53]. This makes it very important to know 
whether this works or not. For some general recent discussions of the scope of holography 
see [51]. Another direction for relating Feynman graphs to string amplitudes, by directly 
identifying the integrand on the moduli space of complex structures of string worlsheets, 
has been advocated in [55J. The ideas emerging from the permutations-strings connection 
may well have interesting overlaps with this programme, which we will leave to future 
investigation. 

If there are indeed new full-fledged gauge-string dualities for theories such as 4 in four 
dimensions, or <p 3 in six dimensions or indeed QED or QED-like theories, there are some 
obvious questions to be addressed. What is the map between physical parameters? In 
the 4> 4 theory, Feynman graphs with v vertices are weighted with g v . We have associated 
them with cylinder worldsheets covering a target space cylinder, with maps of degree Av. 
If the cylinder has area A, the Nambu-Goto action will weight such strings with e~ nA . 
This suggests that the QFT coupling constant g should be identified with e~ A . 

Another obvious question : What is the string coupling of this dual string theory? 
Since all the details of the field theory interaction translate into the group Hi at a bound- 
ary of the cylinder amplitude, interaction in the field theory do not translate directly 
into interactions of the string. A related question is whether there is some modification 
of the Feynman graph counting problem on the field theory side which leads to higher 
genus string worldsheets covering the cylinder. This requires further exploration of the 
sort of string amplitudes that come from various QFT Feynman diagrams. The added 
motivation for these investigations is that the stringy picture leads directly to powerful 
counting results, such as those in Appendix D. 

10.2 Topological Strings and worldsheet methods 

We have argued that there is a link between Feynman graph counting and string ampli- 
tudes by recognizing the combinatorics of worldsheet maps to a cylinder in the Feynman 
graph combinatorics. The string theory in question involves the one that appears in the 
large N expansion of 2dYM with U(N) gauge group. This string theory also has a space- 
time description in terms of topological lattice theory with S n (for all n, hence probably 
better described as S 1 ^) as gauge group. The observables needed to describe Feynman 



51 



graph combinatorics are slightly more general than the ones that appear in U(N) 2dYM. 
They involve boundary observables which are not invariant under the entire S n in the 
lattice TFT description, but only certain specified subgroups H 1 ,H 2 . 

Several worldsheet approaches to the string theory of 2dYM have been proposed 
[HI [9l [10]. It will be very interesting to develop, in these approaches, the boundary 
observables which count the Feynman graphs and their symmetries. The geometrical de- 
scription in terms of gluing world-sheets in section 15.61 should be useful. An interesting 
goal would be to find new worldsheet methods for calculating Feynman graph combina- 
torics, which may be efficient for obtaining asymptotic results. 

10.3 Collect results on orbits 

We have developed a method to think about the combinatorics of Feynman graphs, in- 
volving pairs (So, Si) of combinatoric data. This was discussed from the point of view 
of double cosets, which provide a link to generating functions. From the point of view of 
constructing the Feynman graphs, perhaps the most immediate lesson extracted from the 
study of the pairs (Sq, Si) is that 



A Feynman graph is an orbit of the vertex symmetry group acting on the 
Wick contractions. 



It is a simple exercise to obtain both the vertex symmetry group and the set of Wick 
contractions given a quantum field theory. We have listed this data in the table below for 
the theories that we have considered here. The data is relevant to vacuum graphs - there 
are simple generalizations including external edges. 



Group 


acts on 


Feynman Graph problem 




permutations in [2 2v ] of S 2v 


4 theory 


Si; [S3] 


permutations in [2~2~] of S 3v / 2 


3 theory 


s^fsy x £1,3 [S3] 


permutations in [2 2 ] of S 3v3+ 4 V4 


<p 3 + 4 theory 


S V [S 2 ] 


All permutations in 


Yukawa/QED 


S V [S 2 ] 


All even-cycle permutations in S 2v 


Furry QED 


3 v [Zi] 


permutations in [2 2v ] of S 2v 


Large N expansion of Matrix 4 



The double coset picture allows us to recognize that one can also consider the symmetry 
of the Wick contractions (e.g S^ufSV] in case of <f) A ) acting on the cosets associated with 
the vertices (S^/ S V [S±] in this case). 
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10.4 Galois Group and Feynman Graphs 



We showed that the counting of vacuum graphs in QED/ Yukawa and QED with Furry 
constraint implemented can be mapped to counting of ribbon graphs. In the case of QED 
with Furry constraint, there is an even-only restriction on the vertices. We know, from 
the theory of Dessins d'Enfants [12], that there is an action of the absolute Galois group 
Gal(Q/Q) on ribbon graphs. (For recent applications of Dessins d'Enfants to Hermitian 
Matrix Model correlators, supersymmetric gauge theories and extensive references to the 
Mathematical literature see [T3| [T9| HU [15].) This means that there is an action of the 
Galois group on the QED vacuum diagrams. The list of orders of the vertices is a Galois 
invariant [16], so the Furry theorem restriction continues to allow the Galois action to 
close. 

A related fact is that sums over permutations restricted by conjugacy classes are known 
to be sums over Galois group orbits. 

E E E Kw*) (10.7) 

o-lGTi (T2GT2(T3eT3 

Given the link we have observed between counting ribbon graphs and counting QED 
graphs, we have Galois actions on the QED/Yukawa and Furry QED vacuum graphs. We 
do not know if such is the case for the vacuum graphs of scalar field theories. The answer 
depends on whether sums of the following form 

E E E %^w 3 ) (io.8) 

CTiGHi <7 3 £S n 

for Hi = S v [S±] , H 2 = ifiyfSy are related to Galois theory. We know if Hi = S V [Z^, 
because the connection between QED and ribbon graphs we described, these sums count 
Dessins. For Hi = S V [S^] we do not know an argument to relate to Galois actions. 

11 Summary and Outlook 

We have shown that the counting of Feynman graphs and their symmetry factors in scalar 
field theory, e.g 4 or 3 theories, as well as QED, can be mapped to string amplitudes. 
The string amplitudes are of the type that appear in the string dual of large N Yang 
Mills theory, for which there are several proposed worldsheet constructions. The large N 
2dYM observables are known to be related to S n TFT. The counting problems related to 
Feynman graphs have been expressed as observables in this S n TFT. We used this S n TFT 
data to construct covers of the cylinder. The covers are inrepreted as string world-sheets, 
which are also of cylinder topology. The form of the interactions in the QFT determines 
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the observables at the two ends of the spacetime cylinder, which constrain the windings 
of the string worldsheets to belong to certain subgroups of the permutation group. 

The formulation in terms of string amplitudes is directly related to some classic formu- 
lae on counting of graphs in papers by Read[20]. We have found it useful to think about 
these formulae in terms of an operation of introducing an additional vertex in the middle 
of each edge of the Feynman graph. This separates each edge into a pair of half-edges. 
We label the half-edges with numbers 1, 2, ■ • • , n, and associate a quantity Eo describing 
the vertices of the Feynman graph. The newly-added vertices are described by a quantity 
Ei. This operation is called "cleaning" in the context of ribbon graphs and related Belyi 
maps. In that case two permutations ctq, <J\ play the role of Eo, Ei. In the case at hand, 
E is generically not a permutation, but can be viewed as labeling a coset of permutation 
groups. 

There are groups Hi, H 2 which are symmetry groups of E and Ex respectively. The 
counting of Feynman graphs is equivalent to the counting of elements in the double coset 
of the form 

Hi \ ( Permutation group )/H 2 (H-l) 

For the case of scalar field theories the permutation group is something of the form S n . 
In QED, it is a product group. 

The double coset connects directly to a string amplitude with cylinder target space, 
where Hi, H 2 are associated with the boundaries. The similarities in the current approach 
between Feynman graphs and ribbon graphs, allows us to formulate in group theoretic 
terms, and derive nice formulae for questions such as the total number of types of ribbon 
graphs (summed over genus). Further it allows us to express as a group theory problem 
the counting of the number of ribbon graphs which correspond to the same Feynman 
graph. We may interpret this by saying that large N ribbon graphs (say of matrix (j) 
theory) arise from ordinary Feynman graphs (of 4 theory) by a symmetry breaking of 
S^l/Sy to iS^Z^J. 

The formulae we obtain for QED Feynman graph counting turn out, upon simplifi- 
cation, to be related to the counting of ribbon graphs. We give an explanation of this 
relation by mapping QED graphs to ribbon graphs. The key point is the cyclic nature of 
ribbon graph vertices which map to orientations of fermion loops, in the correspondence 
we describe between ribbon graphs and QED graphs. 

We outline some avenues for extensions of this work. It will be desirable to understand 
in terms of double cosets and string amplitudes and to derive generating functions for 
the refinements of graph counting that occur in QFT. We have made some steps in 
the direction of connected graphs, and a more systematic understanding for the case of 
multiple external legs will be useful. In QFT, it is a familiar fact that the generating 
functional of one-particle irreducible graphs is related to the generateing functional of 
connected graphs by a Legendre transform. Clarifying the implications for double cosets 
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and strings will be desirable. Extending to other quantum field theories is an obvious 
direction. This would give new counting results for the relevant Feynman graphs, and 
could also help address some conceptual issues on the meaning of the string amplitude 
interpretation, with a view to exploring the possibility that the stringy combinatorics 
we have uncovered is merely the tip of an iceberg, the bulk of which is a full-fledged 
QFT-string duality, which does not involve large N. 

With regard to precise information on the asymptotic growth of amplitudes in per- 
turbation theory, the counting sequences of Feynman graphs and their asymptotics are of 
interest. For the case of vacuum graphs in 4 theory, there is an asymptotic result [57] 



e «/4 [fHl! (11 2 ) 

(4!)«(2u)!2 2 *d 1 ' ; 

For general r , with v vertices, we need rv to be even for non- vanishing vacuuum diagrams, 
i.e rv = 2m for some m, and the result [57] is 



-(r 2 -l)/4_N[ 



L m\r\ u v\ 



;ii.3) 



We have given a string interpretation of these sequences. Strings at large quantum num- 
bers often become classical. Can the above asymptotic result, e.g in 4 , be explained by 
an appropriate semi-classical string ? For QED/Yukawa Feynman graphs or for QED, 
with the Furry constraint, we hope that the counting sequences and analytic formulae we 
have described, along with techniques such as those of [57], will allow the determination 
of the asymptotics. 
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A Semi-direct product structure of Feynman graph 
symmetries. 



In this section we will explain that the group of automorphisms of a Feynman graph can 
be realized as the semi-direct product of two subgroups defined shortly. Towards this end, 
it is useful to recall the definition of the semi-direct product. Given two groups G\ and 
G 2 , and a group homomorphism ip : G 2 — > Aut(Gi), the semi-direct product of G\ and 
G2 with respect to ip is denoted G\ Xy, G 2 . As a set G\ xi^ G 2 is the Cartesian product 
Gi x G 2 . Multiplication is defined using ip as 

(gi, hi) * (#2, h 2 ) = (giiph^), hh 2 ) (A.l) 

for all gi, g 2 G Gi and h±, h 2 G G 2 . The identity element e is (e^, ec 2 ) and 

(g,h)- 1 = (1>h-i(9- 1 ),h- 1 ) (A.2) 

The set of group elements (g±, ec 2 ) for a normal subgroup of G\ G 2 isomorphic to G 2 , 
while the set of elements {ec 1 , h) form a subgroup isomorphic to G\. 

To make our discussion concrete, again consider g<j) A theory. Recall from Section H] 
that any graph can be specified, after introducing a new type of vertex (say white when 
the original vertices are colored black) between existing edges and labelling the resulting 
half-edges, by data So associated with the vertices and Si with the edges. For a graph 
with v vertices, So is a collection of v 4-tuples of numbers 

So = ;Qs« ( A -3) 

r=X 

Si is a product of two cycles. The symmetric group S$ v acts by permuting the half edges. 

We can define the Automorphism group of the graph using the data So and Si. It is 
the group of permutations 7 G S^ v which have the property 

7(S ) = S 

7 (Si) 7 ' 1 = Si (A.4) 

By 7(So) we mean the operation which acts on, on each factor of So, as follows 

<ijkl >^< 7 (i) j(k) 7(0 > (A.5) 

The action on Si can also be written as above, or equivalently in terms of conjugation in 
S4,,. In testing the first equality, we treat each angled-bracket as completely symmetric. 

Two types of actions can be automorphisms : vertices (i.e. black dots) can be swapped, 
and propagators (i.e. white dots) can be swapped. Construct the group Ge which acts 



56 



only on the propagators and the group Gy that acts on the vertices. The elements of Gy 
act as a non-trivial permutation on the r index running over the vertices. The elements 
of Ge act trivially on the r index. 

For any given Feynman graph we can argue that Aut(.D) = Ge ><v Gv- Towards this 
end we will prove that Gy acts as an automorphism of Ge- For each element v £ Gy 
there is a permutation 7 £ S v which permutes the vertices 

„(E£>) = 

The elements of Ge leave the vertices fixed 

e(E< r } ) = E« 6 £ G £ 

We have 

u-'eu^) = u-hi^) = u~\{^) = E^ 

Thus, i/ _1 ez/ £ Ge- Defining ijj v {e) = v~ l ei>, it follows that ^ Ul ^ U2 (e) = ip VlU2 (e) which 
shows that we indeed have a homomorphism if) : Gy — > Aut(G ? e). 

To understand why Aut(Z)) = Ge x</>GV, consider the product law for the semi-direct 
product 

(ei,i/i) * (e 2 ,v 2 ) = (ei^(e 2 ), 1/11^2) = (eii^e^ -1 , 1^) 

Acting on the graph D means acting on So and Si as defined above. The left hand 
side applies eiz/ie 2 ^ 2 to D. The right hand side applies ei'0 I/1 (e 2 )z/ii/ 2 = eiVi e 2 z/f x v\y% = 
e\V\€2V2 completing the demonstration. 



B 




Figure 22: For this Feynman graph v — 3. 



An example is in order. For the graph in Figure [221 we have 

S = < 1 2 3 4 >< 5 6 7 8 >< 9 10 11 12 > 

S x = (1 10)(4 11)(3 8)(2 5) (6 9)(7 12) (A.6) 
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From the figure we read off the generators (1 4)(10 11), (2 3) (5 8), (9 12)(6 7) for G E . 
Since these generators commute, Ge is a group of order 8. The group Gy will permute 
the vertices A, B and C of the Feynman graph. There are 6 possible permutations of the 
vertices. To obtain the generators of Gy consider (for example) the permutation which 
swaps B and C, shown in figure |2"31 The relevant permutation is <7bc — (1 2) (3 4) (5 10) (6 
9)(7 12) (8 11). This is indeed an endomorphism of Ge since 

a BC {\ 4)(10 11)<7bo = ( 2 3 )( 5 8 ) <W 2 3) (5 8)a B 1 c = (1 4)(10 11) 

a BC (9 12)(6 7)a B l c = (Q 7) (9 12) 

The same is true of the other elements of Gy. Gy has order 6 and Ge has order 8 so that 
Aut(.D) is order 48. This Feynman graph thus comes with a coefficient (4!) 3 /48 = 288. 




Figure 23: The two Feynman graphs shown are related by swapping vertices B and C. 



Note that there is a notion of Graph automorphism used in graph theory [58]. This 
only includes Gy, since the standard labelling of graphs is to label the vertices and list 
the edges as pairs of numbers associated with the vertices they are incident on. The 
distinction is discussed in [59] . 

B Functions on the double-coset 

In this Appendix we will explain how to build a complete set of functions on S n \ (S n x 
S n ) I (H\ x if 2 ). This Appendix makes use of techniques similar to what was used in [2~T] . 

A basis of functions on S n is given by the matrix elements of irreducible representa- 
tions. Start with an element in the group algebra of (S n x S n ) labelled by representations 
and states in the representations. 
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We can make it invariant under left action of S n and under right action of Hi x H% 
by taking 

= ^im E (« ® «)o£&* ft ® & (b. 2) 

Some basic manipulations lead to 

^ = ^7ii¥T E E E E D hU aa iPi) D ZU aa M n®** 

1 2 o-i,ct 2 g5„ «gs„ ftei^i fee-ff2 

= E EE E ^<(«wt*( a ) fl a(^oy 

*' 111 21 <Ti,<T 2 €Sn aS5„ ft6Hi 
0-1,0-2 Mi,M2 

^Ri^h^^MBh 1 ^ 1 ( 1 H 1 )Bf 1 1 ^\l Hl )B^ 2 ^ 2 (l H2 )Bf^ f12 (1 H2 ) (B.3) 

The label /ii runs over the multiplicity with which the trivial irrep. of H\ appears in the 
decomposition of the S n irrep. R with respect to the subgroup. Thus 



= S RlR2 Y,SniM°- 1 ® C 7 2 )^ 1 ,H 2 A^(^X 1,Ml (l^)^ 1, ' il (l^)< 2 ^(1^)5? ^(1*) 

0,02 

(B.4) 

The sum over o"2 is trivial, so it suffices to consider 

°ln = E a ^hW^H^J^^^i^C 1 *)^^) ( B - 5 ) 

oSS„ 

From this C^^;^ ^ s recons t ruc ted by using Sr^^S^. 

Now define an element in the group algebra of S n labeled by multiplicities (Vi,^) 
of the identity irrep of (Hi,H 2 ) appearing in a decomposition of irrep R of S n . Using 
orthogonality of the branching coefficients we have 

(B.6) 
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So we see that a complete set of functions on the coset 

S n \ (S n X X H 2 ) (B.7) 

can be labelled by the representations R of S n which contain the trivial of H\ and of H 2 . 

We have a complete set of elements in the group algebra of (S n x S n ), which are 
invariant under the left action of S n and the right action of Hi x H2. 

This can be viewed as a derivation of the formula for Feynman graph counting in 
terms of multiplicities given in formulas (17.41) and ( 15.251) . 



B.l QED counting in terms of representation theory 

To count the number of QED vacuum graphs with 2v vertices we need to evaluate 

Fqed(2v) ^ Yl 5 (^7"V) 

' cr£S v [s 2 ] ies 2 v 

v 1 Rr2v a€S v [S 2 ]-f<=S 2v 

1 ' Rr2v a&S v [S 2 \l&S 2v 

= E XrWXrW (B.8) 

' Rh2v a£S v [S 2 ] 

We have used the orthogonality of matrix elements 

E DjkMDg (7- 1 ) = (B.9) 

The product of characters can be be expanded using the Clebsch-Gordan (inner-product) 
multiplicities. C (R, R, A) is the number of times the representation A of appears in 
the tensor product R x R, when this is decomposed in terms of the diagonal S^. Thus 

Fqed(2v) =^yE E C(R,R, A)xa(ct) 

' Rh2v a£S v [S 2 ] 

= J2C(R,R,A)M$ Sv[S2] (B.10) 

R\-2v 

The multiplicity M.^ .„ , is the number of times the identity representation of SUSy ap- 
pears when the representation A of is decomposed into the irreducible representations 
of the subgroup S V [S?\. 

Similar manipulations in the case of Furry QED leads to expressions involving, not 
Clebsch-Gordan multiplicities, but Clebsch-Gordan coefficients, along with the matrix 
elements for ^2 Teven r (g> r in R (g) R. 
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C Feynman graphs with GAP 



We have explained different points of view on the enumeration of Feynman graphs, perhaps 
the most intuitive and useful is that Feynman graph is an orbit of the vertex symmetry 
group acting on the Wick contractions. For the case of 4 theory, we have the wreath 
product S^Sy acting on the conjugacy class of permutations in consisting of 2-cycles, 
which we denoted [2 2v ]. Calculations with this formulation are easy to implement directly 
in the GAP software for group theoretic computations |22j . 

C.l Vacuum graphs of 4 

We illustrate with the sequence 1, 3, 7, 20, 56... of vacuum Feynman graphs in </> 4 theory. 
The count of 3 for v — 2 vertices can be obtained from GAP using the commands. 

gap> C := ConjugacyClass( SymmetricGroup(8), ( 1,2) (3,4) (5,6) ( 7,8) ) 
gap> G := WreathProduct ( SymmetricGroup(4), SymmetricGroup(2) ) 
gap> Length(OrbitsDomain ( G , C ) ) ; 
[3 

This directly implements the formulation of vacuum Feynman graphs as orbits of 
S^SV] on the conjugacy class [2 2v ] described in section HI The numbers 7,20 can also 
be recovered without much trouble with this method. For the purposes of just getting 
the number of vacuum graphs, this method is an overkill. Better methods with cycle 
indices are described in Section [5J We found the use of GAP to be very useful as a tool to 
check in formulating the correct group theoretical formulations of various Feynman graph 
counting problems. The command OrbitsDomain ( G , C ) actually gives a nested list 
of lists of pairings. The list runs over inequivalent Feynman graphs. For each Feynman 
graph, there is a list of different Wick contractions which lead to the given graph. Hence 
these orbits encode not just the number of Feynman graphs, but the Feynman graphs 
themselves. 

C.2 (j) 4 with external edges 

It is a small generalization to consider graphs with external edges. For <ft A theory graphs 
with E external edges we act with the wreath product S.JS4] on the conjugacy class 
[2 2v + t] of Si V+ 2E- The relevant sequence for E = 2 is 1, 2, 7, 23, 85, 340, .... The count of 
7 for v = 2 and E = 2 can be obtained from GAP as follows 

gap> C := ConjugacyClass( SymmetricGroup(lO), ( 1,2) (3,4) (5,6) (7,8) (9,10) ) 
gap> G := WreathProduct ( SymmetricGroup(4), SymmetricGroup(2) ) 
gap> Length(OrbitsDomain ( G , C ) ) ; 
[7 
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Comments made above are again applicable: there are better methods to count these 
graphs which use cycle indices. 

C.3 Symmetry factor 

The symmetry factor of a given Feynman graph can be constructed using GAP. The 
following commands compute the symmetry factor of the diagram in Figure [22j 

gap> HI := WreathProduct( SymmetricGroup(4) , SymmetricGroup(3) ) ; 
gap> Sigl := (1,10) ( 4,11) (3,8) (2,5) (6,9) ( 7,12) ; 
gap> for g in HI do 

if OnPoints( Sigl , g ) := Sigl then 
Count sym := Count sym +1 ; fi ; od ; Count sym ; 

[ 48 

Simple modifications of these commands can construct the group, study its subgroups 
etc. , but these will have no immediate interest for us. 

C.4 Action of S v on a set 

An alternative way to code the graphs is to label the vertices {1, • • ■ , v}. Consider lists 
of 2v unordered pairs. 

{(hji), {i2,32)-~(i2vj2v)} (C.l) 

Put the constraint that 

^6(i,i k ) + 6(i,j k )=4 : (C.2) 

k=l 

for all i from 1, ..v. This imposes the condition that all vertices are 4-valent. 

There is an action of S v on this set. We are interested in the orbits. This can be 
programmed in GAP and gives an alternative method to get the sequence of vacuum 
Feynman graphs in <f) A theory. 

The method of using S v is more efficient in generating the whole set of Feynman 
graphs with GAP, than the S^Sy method. But it does not automatically have the edge 
symmetries, unlike the S^fSV] approach. It only gives the automorphism Gy discussed in 
Appendix [A] 
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D Integer sequences 



In this Appendix we collect the numerical results we have obtained by counting Feynman 
graphs. 

• The number of vacuum graphs in 4 theory, with v vertices, starting with v — 1, is 
given by the sequence 

1, 3, 7, 20, 56, 187, 654, 2705, 12587, 67902, 417065, 2897432, 22382255, 189930004, 
1750561160,.. (D.l) 

This sequence is listed in "The On-Line Encyclopedia of Integer Sequences" (23], 
where it is described as the sequence of 4-regular multi-graphs (loops allowed). This 
sequence is derived using (5.21) or (6.4). 

• The number of connected vacuum graphs in 4 theory with v vertices, starting with 
v = 1, is given by the sequence 

1, 2, 4, 10, 28, 97, 359, 1635, 8296, 48432, 316520, 2305104, 
18428254,... (D.2) 

This sequence was generated using (16. 5p . 

• The number of graphs in 4 theory with E = 2 external legs, with v vertices, starting 
with v — 0, is given by the sequence 

1, 2, 7, 23, 85, 340, 1517, 7489, 41276, 252410, 1706071, 12660012, 
102447112,... (D.3) 

This sequence is generated by the formula (16.81) with E = 2. 

• The number of connected graphs in 4 theory with E = 2 external legs, with v 
vertices, starting with v = 0, is given by the sequence 

1, 1, 3, 10, 39, 174, 853, 4632, 27607, 180148, 1281437, 9896652, 
82610706,... (D.4) 

This sequence is generated by the formula ( 16.91) . These graphs have been tabulated 
in [55]. Comparing the Table II of [55] , we find a match between the number 
of diagrams at first and second order in perturbation theory. At third order in 
perturbation theory we have 10 graphs compared to the 8 graphs listed in [56]. The 
two graphs not listed in [56] are obtained from graph # 5.3 by swapping the position 
of the closed loop and setting sun, and from # 7.2 by swapping the order of the 
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double bubble and the bubble. At fourth order in perturbation theory, there are 30 
graphs listed in [56] compared to our count of 39. However, there are 9 graphs that 
are obtained from graphs #11.2, #11.4, #12.3, #14.2, #14.3, #15.1, #13.3, #15.5 
and #17.2, in much the same way that we described for third order in perturbation 
theory. 

The number of graphs in 4 theory with E — 4 external legs, with v vertices, starting 
with v = 0, is given by the sequence 

3, 10, 44, 190, 889, 4490, 24736, 148722, 976427, 6980529, 54151689, 
453922676,... (D.5) 

This sequence is generated by the formula (16.81) with E — 4. 

The number of vacuum graphs for <p 3 theory, with 2v vertices, starting from v = 1 
2, 8, 31, 140, 722, 4439, 32654, 289519,... (D.6) 
This sequence was generated using (16.121) . 

The number of connected vacuum graphs in <p 3 theory, with 2v vertices, starting 
from v = 1 

2, 5, 17, 71, 388, 2592, 21096, 204638,... (D.7) 
This sequence was generated using (16. 5ft . 

The number of graphs in 3 theory with E = 2 external legs, with 2v vertices 
starting from v = 1, is 

5, 30, 186, 1276, 9828, 86279, 866474, 9924846,... (D.8) 

This sequence was generated using (16.151) with E = 2. 

The number of QED /Yukawa vacuum graphs with 2v vertices, which equals the 
total number of ribbon graphs with 2v edges is given by 

2, 8, 34, 182, 1300, 12634, 153590, 2230979, 37250144,... (D.9) 

This sequence was generated using (18. 27ft . 

The number of connected QED/Yukawa vacuum graphs with 2v vertices, starting 
from v = 1, is given by 

2, 5, 20, 107, 870, 9436, 122832, 1863350, 32019816,... (D.10) 

This sequence was generated using ( 16. 51) . 
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• Vacuum graphs in QED after implementing the constraint due to Furry's theorem. 

1, 4, 12, 57, 321, 2816, 31092, 423947,... (D.ll) 
This sequence was generated using (19.131) . 

• Connected Vacuum graphs in QED after implementing the constraint due to Furry's 
theorem. 

1, 3, 8, 39, 240, 2332, 27196, 382, 802,... (D.12) 
This sequence was generated using (16. 5p . The 1, 3, 8 agree with the results of [43] . 
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